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THERMO-MECHANICAL  AND  THERMAL 
BEHAVIOR  OF  HIGH-TEMPERATURE  STRUCTURAL  MATERIALS 


PREFACE 
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Technical  ceramics  because  of  their  chemical  inertness,  high  melting 
point,  good  wear  resistance,  excellent  mechanical  stability  at  high  tempera¬ 
ture  and  other  unique  properties,  represent  a  class  of  materials  eminently 
suited  for  many  critical  engineering  applications.  Unfortunately,  because 
of  their  brittleness  and  unfavorable  combination  of  pertinent  material 
properties,  technical  ceramics  generally  are  highly  susceptible  to  catas¬ 
trophic  failure  in  non-uniform  thermal  environments,  which  give  rise  to 
thermal  stresses  of  high  magnitude. 

Thermal  stress  failure  analysis  of  structural  materials  represents 
a  multi-disciplinary  problem  which  involves  the  principles  of  heat  trans¬ 
fer,  mechanics  and  materials  engineering.  Over  the  last  few  decades  much 
general  understanding  of  the  nature  of  thermal  stress  failure  of  brittle 
materials  has  been  generated.  However,  due  to  the  multi-disciplinary 
nature  of  the  problem,  the  ability  to  predict  thermal  stress  failure  quan¬ 
titatively  for  design  or  other  purposes  has  lagged  behind  the  progress  made 
in  other  engineering  fields.  The  objective  of  the  present  program  is  to 
improve  the  qualitative  and  quantitative  understanding  of  the  nature  of 
thermal  stress  failure  of  brittle  structural  materials,  including  the  ex¬ 
perimental  as  well  as  theoretical  variables.  In  order  to  achieve  this 
objective,  the  participating  investigators  and  scope  of  the  program  are 
organized  such  that  full  advantage  is  taken  of  the  combined  inputs  from  a 
number  of  engineering  disciplines.  In  a  similar  spirit,  a  number  of  studies 
were  conducted  in  cooperation  with  investigators  at  other  institutions. 

The  effort  of  this  program  consists  of  four  main  themes,  including: 
experimental  thermal  shock  testing  with  supporting  analyses,  measurement 
of  thermophysical  properties  relevant  to  thermal  stress  failure,  the 
analysis  of  mechanisms  of  thermal  stress  failure  and  the  dissemination 
of  information  on  thermal  stresses  in  the  form  of  review  articles,  con¬ 
ferences,  etc. 

A  total  fifty- five  (55)  published,  accepted  or  submitted  technical 
publications  have  resulted  from  this  research  program  since  its  initiation. 
The  status  of  these  publications  is  summarized  in  Table  1. 

Studies  completed  within  the  period  covered  by  this  report  are  pre¬ 
sented  as  individual  chapters  in  the  main  body  of  this  report.  The  title 
of  these  chapters  within  their  main  theme  together  with  a  brief  comment 
are  as  follows: 
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A.  Measurement  of  Thermophysical  Properties  Relevant  to  Thermal  Stress 
Failure. 

Chapter  I:  Y.  Chen,  M.  M.  Abraham,  L.  D.  Bentsen  and  D.  P.  H.  Hassel- 
man,  "Effect  of  Ni-Alloyinq  on  the  Thermal  Diffusivity/Con- 
ductivity  of  MgO  Single  Crystals" 

This  study  showed  that  the  thermal  diffusivity/ conductivity  of  single 
crystals  of  magnesium  oxide  can  be  affected  stronqly  by  the  alloying  with  Ni. 
The  relative  effect  was  found  to  be  a  function  of  whether  the  Ni  existed  as 
a  solid-solution  or  in  the  form  of  precipitated  colloids. 

Chapter  II:  L.  D.  Bentsen,  D.  P.  H.  Hasselman  and  R.  Ruh,  "Effect  of 
Hot-Pressing  Temperature  on  the  Thermal  Diffusivity/ Conductivity  of  SiC-AIN 
Composites" 

The  results  of  this  study  show  that  the  heat  conduction  of  mixtures  of 
materials  depend  critically  on  whether  such  mixtures  exist  in  the  form  of 
solid  solutions  or  as  distinct  composites.  Specifically,  it  was  found  that 
the  thermal  conductivity/diffusivity  of  SiC-AIN  mixtures  can  be  controlled 
by  varying  the  hot-pressing  temperature  which  in  turn  controls  the  degree  of 
inter-diffusion  between  the  Si C  and  AIN  components. 

Chapter  III:  M.  A.  Bucknam,  L.  D.  Bentsen  and  D.  P.  H.  Hasselman, 

"The  Measurement  of  the  Thermal  Conductivity  of  Refractories  by  the  Laser- 
Flash  Method" 

This  feasibility  study  showed  that  the  laser- flash  method  which  uses 
dime-sized  specimens  can  be  used  successfully  to  measure  the  temperature 
dependence  of  the  thermal  diffusivity  and  conductivity  of  materials  with 
relatively  coarse  micro-structures,  as  long  as  care  is  exercised  in  se¬ 
lecting  a  specimen  with  a  thermal  diffusivity  as  close  as  possible  to  the 
mean  of  the  value  for  a  batch  of  specimens.  Data  obtained  for  a  series  of 
refractories  selected  as  appropriate  materials  for  this  study  showed  ex¬ 
cellent  agreement  with  the  thermal  conductivity  obtained  by  two  other  inde¬ 
pendent  methods. 

B.  Analyses  of  Role  of  Material  Properties  in  Thermal  Stress  Resis¬ 
tance  of  Brittle  Structural  Ceramics 

Chapters  IV,  V,  VI,  VII  and  VIII: 

J.  R.  Thomas,  J.  I.  Frankel  and  D.  P.  H.  Hasselman,  "Effect  of  Inter¬ 
face  Reflections  and  Angle  of  Incidence  of  Radiation  on  Thermal  Stresses 
in  Semi-Transparent  Materials," 

J.  P.  Singh,  J.  R,  Thomas,  Jr.,  and  D.  P.  H.  Hasselman,  "Role  of  Ab¬ 
sorption  Coefficient  in  the  Frequency  Dependence  of  the  Thermal  Stresses 
in  a  Partially-Absorbing  Plate  Subjected  to  Cyclic  Thermal  Radiation," 

J.  P.  Singh,  J.  I.  Frankel,  J.  R.  Thomas,  Jr.,  and  D.  P.  H.  Hasselman, 
"Thermal  Stresses  in  Partially  Absorbing  Flat  Plate  Due  to  Sudden  Inter¬ 
ruption  of  Steady-State  Asynmetric  Radiation,  I:  Convective  Cooling  at 
Rear  Surface," 
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J.  I.  Frankel ,  J.  P.  Singh,  J.  R.  Thomas,  Jr.,  and  D.  P.  H.  Hasselman, 
"Thermal  Stresses  in  Partially  Absorbing  Flat  Plate  Due  to  Sudden  Inter¬ 
ruption  of  Steady- State  Asymmetric  Radiation,  II:  Convective  Cooling  at 
Front  Surface," 

J.  P.  Singh,  J.  R.  Thomas,  Jr.,  and  D.  P.  H.  Hasselman,  "Thermal  Stresses 
in  a  Partially  Absorbing  Flat  Plate  Asymmetrically  Heated  by  Cyclic  Thermal 
Radiation  and  Cooled  by  Convection:  Addendum." 

These  chapters  represent  a  continuation  of  the  analysis  of  the  role  of 
the  absorption  coefficient  in  the  thermal  stress  failure  of  partially  absorb¬ 
ing  ceramic  materials  subjected  to  thermal  radiation  for  a  variety  of  boundary 
conditions.  In  general,  the  analytical  results  show  that  the  absorption  coef¬ 
ficient  plays  a  significant  role  in  governing  the  magnitude  of  thermal  stress, 
coupled  with  the  frequency,  angle  of  incidence,  multiple  reflections  of  the 
incoming  radiation  and  nature  of  convective  cooling.  This  phase  of  the  pro¬ 
gram  has  been  completed. 

Chapter  IX:  W.  Chang,  C.  E.  Knight,  D.  P.  H.  Hasselman  and  R.  G.  Mitchi- 
ner,  "Analysis  of  Thermal  Stress  Failure  of  Segmented  Refrac¬ 
tory  Structures" 

This  analysis  concentrated  on  the  role  of  the  material  properties  and 
geometry  in  a  unique  failure  mode  encountered  in  segmented  high- temperature 
structures,  which  consists  of  crack  formation  perpendicular  to  the  direction 
of  heat  flow  in  contrast  to  the  more  common  mode  of  crack  propagation  parallel 
to  the  direction  of  heat  flow.  For  other  specific  boundary  conditions  chosen, 
it  was  found  that  high  values  of  thermal  diffusivity  are  desirable  if  the 
magnitude  of  thermal  stresses  are  to  be  minimized.  Furthermore,  the  magnitude 
of  the  dimensions  of  the  individual  segments  also  are  critical  in  qoverning 
the  value  of  thermal  stress. 

C.  Reviews 

Chapter  X:  D.  P.  H.  Hasselman  and  L.  D.  Bentsen,  "Effect  of  Microstruc- 
tural  and  compositional  Heterogeneity  on  the  Conduction  of 
Heat  in  Structural  Materials  for  High-Temperature  Use" 

This  review  article  represents  a  suronary  of  earlier  results  obtained  for 
the  compositional  and  microstructural  effects  on  the  conduction  of  heat  in 
structural  ceramic  materials. 
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Steady  State  Heat  Flow,"  J.  Amer.  Ceram.  Soc.,  62,  432  (1979). 
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EFFECT  OF  Ni -ALLOYING  ON  THE  THERMAL 


DIFFUSIVITY/CONDUCTIVITY  OF  MgO  SINGLE  CRYSTALS 
Y.  Chen,  M.  M.  Abraham,  L.  D.  Bentsen,  D.  P.  H.  Hasselman 

The  thermal  diffusivity/conductivity  of  Ni-doped  MgO  single  crystals 

was  measured  as  a  function  of  temperature  by  the  laser-flash  method.  The 

2+ 

Ni  in  the  form  of  Ni  substitutional  solid  solution  has  a  much  greater 
effect  on  the  thermal  diffusivity/conductivity  than  precipitated  Ni  colloids. 
These  results  are  discussed  in  terms  of  phonon  and  photon  scattering  effects. 


The  conduction  of  heat  through  solids  can  be  affected  significantly 

by  the  presence  of  other  elements  and  is  particularly  sensitive  to  their 

size  and  distribution.  Isolated  point  defects  reduce  the  thermal  con- 

1-3 

ductivity  significantly  due  to  increased  phonon  scattering  .  Larger 

colloids  produce  a  similar  effect  but  the  reduction  in  thermal  conductivity 

is  not  as  great  as  for  point  defects.  Still  larger  particles  require  that 

composite  theory  be  used  to  describe  the  overall  thermal  conductivity  which 

4 

may  be  higher  or  lower  than  that  of  the  host  matrix  . 

Studies  of  the  effect  of  metallic  colloids  on  thermal  conductivity 

have  concentrated  mostly  on  the  halides  due  to  the  ease  by  which  they  can 
5  6 

be  synthesized  ’  .  For  refractory  oxides,  Li,  Na  and  k  atoms  have  been 

O  Q 

introduced  by  ion  implantation  ,  however,  these  elements  were  limited  to 
the  surface  regions.  Such  materials  were  therefore  not  suitable  for  measurement 
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ment  of  the  effect  of  alloying  elements  on  the  heat  conduction  behavior 
of  refractory  oxides.  Recently  a  method  which  produces  a  more  uniform 
distribution  of  the  alloying  element  was  developed  and  applied  to  the 
alloying  of  magnesium  oxide  with  NiJ0,ll,12  The  aim  of  the  present  study 
was  to  determine  the  effect  of  the  distribution  of  the  Ni  atoms  on  the 
thermal  diffusivity/conductivity  of  MgO  single  crystals. 

1 3 

The  Ni-doped  MgO  crystals  were  grown  by  the  arc-fusion  technique 

from  high-purity  MgO  powder  with  5  vit.%  NiO.  The  as-grown  crystals  ex- 

2+ 

hibited  a  slightly  greenish  color  associated  with  Ni  in  substitutional 
14  15 

sites  *  .  Atomic  emission  spectroscopy  showed  the  nickel  concentration 

to  be  ^  0.4  wt.%.  Transmission  electron  microscopy  and  optical  absorption 

O  .  I  £ 

measurements  confirmed  that  the  nickel  consisted  primarily  of  Ni  ions  . 
Precipitation  of  the  nickel  was  achieved  by  heating  the  Ni-doped  crystals 
at  210UK  in  a  reducing  atmosphere  of  magnesium  vapor  at  ^  6  atmospheres 
for  periods  of  2  or  20  hours.  Following  reduction  the  Ni-doped  crystals 
appeared  to  be  black. 

Optical  absorption  spectra  at  77K  before  and  after  reduction  of  a 

typical  specimen  is  shown  in  Fig.  1.  The  absorption  spectrum  before 

2+ 

reduction  results  from  the  Ni  ions.  The  large  absorption  coefficient 

of  'v  600  cnf^  in  the  visible  region  is  indicative  of  the  high  opacity 

after  reduction.  The  broad  band  of  high  absorption  coefficient  near  2  eV 

is  attributed  to  the  formation  of  the  Ni  colloids. 

Figure  2a  shows  the  EPR  spectrum  at  80K  of  an  "as-grown"  MgO:Ni  crystal 

with  the  magnetic  field  ranging  from  0  to  500  Gauss  oriented  along  the  <100> 

2+ 

axis.  The  single  isotropic  line  is  due  to  the  Ni  species.  Figure  2b 
shows  the  corresponding  spectrum  of  a  sample  with  the  same  volume,  but 


2+ 

reduced  for  2  hours.  The  Ni  signal  is  markedly  reduced  in  intensity 
and  a  large  broad  background  signal  has  appeared.  The  latter  is  attributed 
to  the  metallic  nickel  colloids  which  grew  at  the  expense  of  the  substitutional 
nickel  ions.  The  sample  reduced  for  20  hours  exhibited  the  broad  spectrum 

24- 

due  to  the  colloidal  nickel  without  evidence  for  the  presence  of  Ni  ions. 

Figure  3  shows  the  data  for  the  thermal  diffusivity  measured  by  the 
laser- flash  method^  for  a  pure  MgO  single  crystal  and  the  single  crystals 
doped  with  Ni  and  heat-treated  as  described  above.  The  pure  crystal  generally 
exhibits  the  highest  values  for  the  thermal  diffusivity,  whereas  the  crystal 
with  nickel  in  solid  solution  exhibits  the  lowest  values.  Precipitation  of 
the  Ni  leads  to  an  increase  of  the  thermal  diffusivity  with  the  heat-treatment 
for  20  hours  resulting  in  values  nearly  identical  to  that  of  the  undoped 
crystal.  Similar  trends  are  apparent  in  the  data  for  the  thermal  con¬ 
ductivity,  shown  in  Fig.  4  as  a  function  of  1/T,  calculated  from  the 

experimental  data  for  the  thermal  diffusivity,  the  density  of  the  crystals 

1 8 

and  handbook  data  for  the  specific  heat  of  magnesium  oxide  and  nickel 
The  linear  dependence  of  the  thermal  conductivity  on  1/T  for  the  Ni-doped 
crystals  indicates  that  heat  transfer  through  the  crystals  occurs  primarily 
by  phonon  transport.  The  positive  deviation  at  higher  temperatures  for  the 
undoped  crystal  suggests  that  heat  transfer  by  radiation  makes  a  significant 
contribution  to  the  total  conductivity. 

In  order  to  provide  an  explanation  for  the  present  observations  it 
may  be  noted  that  the  phonon  conductivity  (K^)  can  be  expressed  as: 

KL  «  1/3  /  c(u)vl(tt)dM  0) 

where  c  is  the  specific  heat  per  unit  volume,  v  is  the  average  velocity 


of  the  phonons,  l  is  the  phonon  mean  free  path,  and  u>  is  the  phonon  fre¬ 
quency.  The  total  mean  free  path  is  composed  of  an  intrinsic  term  )  result¬ 
ing  from  anharmonic  processes  and  a  term  (t^)  resulting  from  various  defects 
wi  th : 


Ml  =  1/^  +  1/Ad  (2) 

For  the  present  samples  the  defect  mean  free  path  is  governed  by  the 
Ni1-'  point  defects  and  the  nickel  colloids.  For  the  same  number  of  Ni 
atoms,  the  probability  of  phonon  scattering  with  corresponding  value  of  l/td 
is  much  higher  for  Ni  atoms  in  solid  solution  than  the  corresponding  value 


,.2+ 


for  extended  Ni  colloids.  For  this  reason,  the  value  for  the  thermal  con- 

„2+ 


ductivity  of  MgO  with  Ni  atoms  is  expected  to  be  well  below  the  corresponding 
value  for  the  MgO  with  Ni  colloids,  in  general  agreement  with  the  present  data. 

It  may  be  noted  that  the  pure  MgO  crystal  has  a  thermal  conductivity 
somewhat  below  the  value  for  the  MgO  with  Ni  colloids.  Possibly  this  observa¬ 
tion  is  attributable  to  a  small  contribution  of  heat  transfer  by  electron 
transport  within  the  Ni  colloids.  However,  this  contribution  is  sufficiently 
small  that  the  overall  conductivity  is  controlled  by  phonon  transport  within  the 
MgO  as  indicated  by  the  1/T  dependence.  This  latter  observation  also  Indicates 


that  the  radiative  component  of  the  total  conductivity  is  also  very  small. 

.19. 


The  thermal  conductivity  due  to  radiation  can  be  expressed  as  either 


(r  55  T  an2r3j 


or 


KR  =  (2?¥  )  4on2l3d 


(3) 

(4) 


where  a,  n,  T,  l ,  d  and  e  are  the  Stefan-Boltzman  constant.  Index  of  re¬ 
fraction,  temperature,  photon  mean  free  path,  sample  thickness,  and  the 


emissivity  of  the  specimen  surfaces,  respectively.  Equation  3  is  applicable 
when  t,  which  is  equal  to  the  reciprocal  of  the  absorption  coefficient,  is 
small  compared  to  the  thickness  of  the  sample.  Such  is  the  case  for  the  Ni 
colloid  containing  crystal  with  an  absorption  coefficient  of  =  600  cm’^  in 
which  less  than  1$  of  the  total  conductivity  may  be  attributed  to  photon 
transfer.  Equation  4  may  be  used  when  i  »  d;  however,  since  the  laser-flash 
method  employs  very  small  samples  (d  8  2  mm)  the  radiative  conductivity  is 
also  small  even  for  the  pure  MgO  crystal ,  except  at  very  high  temperatures, 
in  agreement  with  the  present  data. 

The  data  and  discussion  presented  thus  confirm  the  critical  role  of  the 
distribution  of  alloying  elements  on  the  conduction  of  heat  through  soltds, 
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ABSTRACT 

Data  are  presented  for  the  thermal  dif fusivity/conductivity 
of  hot-pressed  mixtures  of  SiC  and  AIN.  The  results  indicate 
that  hot-pressing  at  higher  temperatures,  which  permits  solid-so¬ 
lution  formation,  results  in  significantly  lower  thermal  diffu- 
sivity/conductivity  than  obtained  by  hot-pressing  at  lower  temp¬ 
eratures,  at  which  the  SiC-AIN  exists  as  discrete  phases. 
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The  thermal  conductivity  of  solids  can  be  affected  signifi¬ 
cantly  by  the  presence  of  foreign  elements  and  is  particularly 
sensitive  to  their  size  and  distribution.  Due  to  increased  pho¬ 
non  scattering,  isolated  point  defects  can  significantly  reduce 
the  thermal  conductivity  (1-3).  In  contrast,  particles  which  are 
large  enough  to  retain  their  individual  bulk  properties  may  cause 
an  increase  or  decrease  in  overall  thermal  conductivity,  depend¬ 
ing  on  their  thermal  conductivity  relative  to  the  host  matrix. 

In  this  case,  thermal  conductivity  can  be  calculated  by  means  of 
composite  theory  (4).  Th.  important  role  of  the  distribution  of 
a  foreign  element  on  heat  conduction  was  demonstrated  recently  by 
experimental  data  (5)  for  MgO  single  crystals  alloyed  with  Ni, 
which  existed  in  either  solid  solution  or  in  the  form  of  precipi¬ 
tates  . 

Because  of  its  high  sensitivity  on  the  distribution  of  for¬ 
eign  elements,  the  thermal  conductivity  of  solids  can  be  cont¬ 
rolled  by  modifying  the  processing  conditions  which  affect  the 
distribution  of  the  alloying  elements.  The  feasibility  of  this 
approach  is  demonstrated  by  experimental  data  for  the  effect  of 
aluminum  nitride  additions  on  the  thermal  diffusivity  of  silicon 
carbide . 

The  specimens  used  for  this  study  were  prepared  for  an  ear¬ 
lier  investigation  of  the  crystallographic,  compositional  and 
microstructural  analysis  of  SiC-AIN  solid  solutions  and  two-phase 
mixtures  (6).  These  samples  were  hot-pressed  from  mixtures  of 
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SiC  and  AIN  powders  at  temperatures  ranging  from  1700  to  2200° C. 
The  specimens  hot-pressed  at  the  lower  temperatures  consisted  of 
mixtures  of  SiC  and  AIN  grains  with  only  slight  inter-diffusion 
between  the  SiC  and  AIN.  These  specimens  prepared  at  the  highest 
temperatures  consisted  of  solid-solutions  due  to  inter-diffusion 
between  the  two  components.  Although  some  of  these  solid-solu¬ 
tions  consisted  of  a  single  polymorph,  the  distribution  of  the 
individual  atomic  species  was  found  to  be  highly  non-uniform. 

For  details  the  reader  is  referred  to  the  original  study  (6). 

The  thermal  diffusivity  of  the  samples  was  measured  by  the 
flash  method  (7)  using  a  glass-Nd  laser  as  the  flash  source,  with 
an  estimated  accuracy  of  ^  3%.  Figure  1  shows  the  experimental 
data  for  the  thermal  diffusivity  at  room  temperature  as  a  func¬ 
tion  of  AIN  content  and  a  range  of  hot-pressing  temperatures.  It 
is  evident  that  the  AlN-SiC  mixtures  hot-pressed  at  the  higher 
temperatures  exhibit  a  thermal  diffusivity  significantly  lower 
than  those  hot-pressed  at  the  lower  temperatures.  The  differ- 
ences  in  the  thermal  diffusivity  for  the  two  AIN  samples  without 
SiC  at  least  in  part  are  due  in  part  to  differences  in  density 
which  was  95.3%  for  the  sample  hot-pressed  at  1700°C  and  94.4% 
for  the  sample  prepared  at  2100° C.  Undoubtedly  other  variables 
may  have  contributed  to  these  differences  as  well.  For  the  SiC 
samples  without  AIN  the  values  for  the  thermal  diffusivity  are 
appreciably  lower  than  the  value  near  0.8  cm  s  observed  in  this 
laboratory  (8)  for  other  samples  of  SiC  from  other  sources.  The 
lower  values  of  this  study  are  thought  to  be  due  to  the  rela- 


4 


tively  high  boron  content  of  the  two  starting  powders  used  for 
their  preparation.  This  hypothesis  is  based  on  preliminary  data 
obtained  during  another  study  (8)  which  showed  that  a  15%  porous 
SiC  without  boron  exhibited  a  thermal  diffusivity  higher  than 
fully  dense  SiC  with  boron. 

Figure  2  shows  the  values  of  the  thermal  conductivity  calcu¬ 
lated  from  the  data  for  the  thermal  diffusivity  given  in  Fig.  1 
by  multiplication  with  the  product  of  density  and  specific  heat. 
Also  included  in  Fig.  2  are  the  data  for  the  thermal  conductivity 
reported  by  Rafaniello  et  al  (9)  for  SiC-AIN  solid  solutions  made 
by  carbothermal  reduction  of  silica  and  alumina  followed  by  hot- 
pressing  at  1950  to  2200 °C.  These  latter  data  show  quite  good 
agreement  with  the  samples  of  the  present  study  hot-pressed  at 
the  two  highest  temperatures  of  2100  and  2200 °C.  Included  in 
Fig.  2  is  the  expected  dependence  of  the  thermal  conductivity 
calculated  from  composite  theory  on  the  assumption  that  the  A1N- 
SiC  consists  of  a  mixture  of  AIN  and  SiC  without  any  inter-diffu¬ 
sion  or  interfacial  resistance  to  the  flow  of  heat.  For  clarity, 
the  composite  thermal  conductivity  was  calculated  from  the  data 
for  SiC  and  AIN  hot-pressed  at  2100 °C.  Clearly,  the  composite 
behavior  will  vary  with  the  choice  of  the  particular  end-point 
data.  Comparison  of  the  calculated  data  for  composite  behavior 
with  the  experimental  data,  especially  those  for  the  highest 
hot-pressing  temperatures,  clearly  shows  that  the  formation  of  a 
solid  solution  results  in  a  dramatic  decrease  of  the  thermal  con¬ 
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The  thermal  conductivity  of  the  AlN-SiC  mixtures  hot-pressed 
at  temperatures  ranging  from  1850  to  2000°C  is  significantly 
higher  than  that  of  the  samples  prepared  at  2100  or  2200°C.  This 
difference  indicates  the  absence  of  complete  solid  solution. 
Presumably,  at  these  lower  temperatures  the  rate  of  inter-diffu¬ 
sion  is  too  low  for  solid-solutions  to  form  so  that  in  these  sam¬ 
ples  the  SiC  and  AIN  have  retained  their  original  form.  It  is  of 
interest  to  note,  however,  that  the  data  for  the  thermal  conduc¬ 
tivity  of  these  samples  falls  below  the  calculated  behavior  if 
the  SiC-AIN  mixtures  were  true  composites.  This  suggests  the 
existence  of  either  a  barrier  to  heat  flow  at  the  SiC-AIN  inter¬ 
face  due  to  less  than  perfect  atomic  bonding,  or  alternatively, 
and  perhaps  more  likely,  solid  solution  formation  occurred  at  the 
interfacial  zone  over  a  thickness  sufficient  to  affect  the  ther¬ 
mal  conductivity  significantly.  In  fact,  such  a  thin  zone  of 
inter-diffusion  was  detected  by  XRD  in  TEM. 

At  higher  temperatures,  where  phonon  scattering  is  greater, 
the  effect  of  the  distribution  of  other  elements  on  thermal  con¬ 
ductivity  becomes  less  significant  (2).  This  behavior  is  exhib¬ 
ited  by  the  present  samples  as  shown  in  Fig.  3  which  compares  the 
temperature  dependence  of  the  thermal  diffusivity  of  two  samples 
containing  50%  AIN.  These  samples  were  hot-pressed  at  1900  and 
2100°C  and  constitute  a  two-phase  mixture  and  a  solid  solution 
respectively.  •  The  relative  difference  in  thermal  diffusivity 
decreases  with  increasing  temperature.  The  same  effect  is  also 
shown  in  Fig.  4  for  two  samples  with  24.8  and  34%  AIN  hot-pressed 
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at  1950  and  2100 °C,  respectively.  Nevertheless,  even  at  the 
higher  temperatures,  significant  differences  in  the  thermal  dif- 
fusivity  of  the  mixtures  and  solid-solutions  still  exist. 

The  results  of  the  present  study  clearly  indicate  that  by 
modifying  processing  conditions  (temperature  in  this  case),  a 
significant  measure  of  control  can  be  exerted  over  the  thermal 
transport  behavior  of  composites  that  form  solid  solutions. 
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Fig.  2.  Calculated  values  of  the  thermal 

conductivity  of  SiC-AIN  composites 
hot-pressed  over  a  range  of  tempera 
ture,  compared  with  literature  data 


TEMPERATURE  (°C) 


Temperature  dependence  of  the  thermal 
diffusivity  of  SiC-AlN  mixtures  with 
50  wt.%  AIN,  hot-pressed  at  1900  and 
2100°C. 
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ABSTRACT 

A  study  was  conducted  to  determine  the  feasibility  of  the 
laser-flash  diffusivity  method  for  determining  the  thermal  con¬ 
ductivity  of  refractories.  Considerable  specimen-to-specimen 
variation  in  the  values  of  thermal  diffusivity  at  room  tempera¬ 
ture  was  found  for  a  set  of  specimens  cut  from  eleven  different 
refractories.  The  specimen  with  the  value  of  thermal  diffusivity 
closest  to  the  mean  was  taken  to  be  representative  of  the  bulk 
material  and  was  used  for  measurements  at  higher  temperatures. 

For  the  dense  refractories  the  thermal  conductivity  at  elevated 
temperatures  obtained  by  the  laser-flash  method  agreed  quite  well 
with  the  corresponding  data  obtained  by  the  guarded-plate  and 
hot-wire  methods.  For  a  relatively  porous  insulating  refractory 
the  thermal  conductivity  obtained  by  the  laser-flash  method  was 
found  to  exceed  the  data  obtained  by  the  guarded-plate  method. 
This  observation  was  thought  to  be  due  to  penetration  of  the  las¬ 
er-flash  into  the  specimen  interior.  The  thermal  conductivity 
near  room  temperature  obtained  by  the  laser- flash  method  exceeded 
the  data  obtained  by  the  guarded-plate  or  hot-wire  methods. 
Further  research  is  needed  in  order  to  establish  the  reasons  for 
this  effect.  In  general,  the  results  of  this  study  indicate 
that,  provided  care  is  taken  in  selecting  the  appropriate  test 
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specimen,  the  laser-flash  method  represents  a  valuable,  conve- 
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1.  INTRODUCTION 

Thermal  conductivity  plays  a  critical  role  in  the  perfor¬ 
mance  of  refractory  materials  in  structures  such  as  high  tempera¬ 
ture  process  vessels,  furnaces,  kilns,  retorts,  and  recuperators 
to  name  just  a  few.  Low  values  of  thermal  conductivity  are 
required  if  heat  losses  are  to  be  kept  to  a  minimum.  On  the 
other  hand,  heat  transfer  through  brickwork  (desirable  in  certain 
heating  operations)  and  control  of  hot-face  temperatures  is 
obtained  most  easily  by  using  refractories  of  higher  thermal  con¬ 
ductivity  (1).  This  latter  characteristic  is  also  desirable,  if 
not  critical,  in  avoiding  the  incidence  of  failure  by  thermal 
shock.  For  these  reasons,  reliable  thermal  conductivity  data  are 
essential  for  the  selection  of  a  refractory  which  will  perform 
optimally  in  a  given  application. 

The  measurement  of  the  thermal  conductivity  of  refractories 
at  this  time  relies  primarily  on  the  guarded-plate  and  hot-wire 
methods  (2,3,4).  An  increasingly  popular  method  for  measuring 
the  conduction  of  heat  through  solids  is  the  laser- flash  method, 
pioneered  by  Parker,  et.  al.  (5).  In  general,  this  method  con¬ 
sists  of  subjecting  one  side  of  a  small,  thin,  disc-shaped  speci¬ 
men  to  a  uniform  laser  pulse  (or  other  suitable  instantaneous 
heat  source).  The  heat  absorbed  at  the  irradiated  surface  will 
diffuse  through  the  specimen.  The  transient  temperature  response 
at  the  opposite  face,  recorded  by  an  appropriate  method,  allows 
calculation  of  the  thermal  diffusivity.  Thermal  conductivity  can 
be  calculated  by  multiplying  the  thermal  diffusivity  with  inde¬ 
pendently  obtained  values  of  density  and  specific  heat. 
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The  specimen  required  for  the  laser-flash  method  is  of 
simple  geometry,  with  dimensions  of  the  order  of  2mm  in  thickness 
by  10mm  in  diameter  or  diagonal,  for  specimens  of  circular  or 
square  geometry.  One  essential  requirement  for  the  measurement 
of  thermal  conductivity  (or  any  property)  by  any  method  is  that 
the  test  specimen  used  for  the  measurement  must  be  representative 
of  the  bulk  of  the  material  under  investigation.  For  specimens 
of  the  dimensions  given  above  for  the  laser-flash  method,  this 
requirement  is  certainly  met  if  the  test  material  is  homogeneous 
and  consists  of  a  single  solid  phase  with  grain  and  pore  sizes  of 
a  few  microns. 

Refractories  are  usually  highly  heterogeneous,  both  chemi¬ 
cally  as  well  as  structurally.  Furthermore,  refractories  possess 
a  relatively  coarse  microstructure  with  a  wide  distribution  of 
pore  and  grain  sizes,  perhaps  as  large  as  a  few  millimeters. 

This  latter  dimension  is  of  the  same  order  as  the  thickness  of 
the  laser  flash  specimens.  Thus,  the  data  obtained  by  the  laser- 
flash  method  may  depend  critically  on  the  location  in  a  refrac¬ 
tory  brick  from  which  the  specimen  was  cut.  Considerable  speci¬ 
men  to  specimen  variation  is  expected  and  may  pose  a  limitation 
on  the  use  of  the  laser-flash  method  for  determining  the  thermal 
dif fusivity/conductivity  of  refractories.  The  purpose  of  this 
study  was  to  resolve  this  question  by  comparing  data  for  the 
thermal  conductivity  of  a  number  of  refractories  measured  by  the 
laser-flash  method  with  corresponding  data  obtained  by  the  guard¬ 
ed-plate  or  hot-wire  method. 
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2 .  EXPERIMENTAL 

2.1  MATERIALS  AND  SAMPLE  PREPARATION 

Table  1  lists  the  eleven  refractory  compositions  selected 
for  this  study.  The  listed  value  of  density  represents  the  mean 
of  the  range  of  density  reported  by  the  manufacturers. 

In  preparing  the  disc-shaped  specimens  required  for  the  las¬ 
er-flash  method,  cores  approximately  12.5mm  in  diameter  by 
approximately  50mm  long  were  cut  from  the  larger  samples  by  means 
of  a  diamond  core  drill.  Because  of  possible  preferred  alignment 
of  the  grains  and/or  pores,  the  heat  conduction  of  refractories 
could  exhibit  anisotropic  behavior.  In  order  to  assure  a  valid 
comparison  of  the  data,  the  cores  from  the  refractory  samples 
provided  by  manufacturer  A  were  drilled  in  the  direction  of  heat 
flow  used  in  the  guarded  plate  method.  Since  the  hot-wire  method 
involves  radial  heat  flow,  the  cores  cut  from  the  samples  pro¬ 
vided  by  manufacturer  B  were  cut  in  two  mutually  perpendicular 
directions,  each  in  turn  perpendicular  to  the  length  of  the  wire 
in  the  hot-wire  method.  The  cores  were  then  sliced  into  a  number 
of  circular  discs  approximately  2mm  in  thickness  with  a  low-speed 
diamond  banch-saw.  The  disc  samples  were  then  cleaned  in  acetone 
to  remove  the  saw  lubricant  and  cutting  debris,  followed  by  dry¬ 
ing  for  at  least  24  hours. 
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2.2  EXPERIMENTAL  PROCEDURE 

2.2.1  Laser-Flash  Method 

Figure  1  schematically  illustrates  the  principle  components 
of  the  experimental  equipment  used  for  this  study.  The  laser- 
flash  source  consists  of  a  Nd-glass  laser*  with  a  laser  rod  hav¬ 
ing  a  diameter  of  approximately  12.5mm.  The  wavelength  of  the 
laser  pulse  was  1.06  vim  with  a  pulse-width  of  about  800  usee.  The 
energy  per  pulse  was  a  maximum  of  approximately  70  Joules  at  the 
5kV  peak  voltage  of  the  laser  power  supply.  The  latter  voltage 
was  kept  at  approximately  4.1kV  in  order  to  extend  the  life  of 
the  flash-lamp.  At  this  voltage,  the  intensity  of  the  beam  was 
highly  uniform,  as  verified  by  the  flash  imprint  on  photographic 
film. 

For  measurement  of  the  thermal  diffusivity  at  elevated  temp¬ 
eratures  the  sample  rested  upon  a  graphite  specimen  holder  which 
consisted  of  a  graphite  ring  6mm  thick  with  an  outer  and  inner 
diameter  of  24  and  6mm,  respectively.  The  specimen  was  centered 
on  a  circular  knife-edge  machined  at  the  edge  of  the  center  hole. 
The  specimen  holder  and  sample  were  held  within  a  carbon  resis¬ 
tance  furnace**  which  was  capable  of  reaching  temperatures  as 
high  as  2500°C  in  an  inert  atmosphere.  Nitrogen  gas  was  used 
throughout  this  study.  The  furnace  was  oriented  vertically  with 
the  specimen  lying  horizontally  within  the  holder.  The  flash 
from  the  laser  entered  the  top  of  the  furnace  via  a  first-surface 


*  Model  K-l,Compulaser,  Canoga  Park,  California 

**Model  1000A,  Astro  Industries,  Santa  Barbara,  California. 
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gold  mirror  through  a  fused  quartz  window.  Because  of  a  slight 
divergence  of  the  laser  beam,  specimens  with  the  same  nominal 
diameter  as  the  laser  rod  were  irradiated  uniformly  over  their 
entire  surface,  also  verified  by  the  flash  imprint  on  photo¬ 
graphic  film. 

The  transient  temperatures  of  the  lower  face  of  the  specimen 
induced  by  the  laser  pulse  at  the  top  surface  were  monitored  by 
remote  optical  sensors.  An  InSb  inf rared-detector*  was  used  for 
specimen  temperatures  ranging  from  about  300° C  to  600° C.  For 
temperatures  in  excess  of  600 °C  the  transient  temperature  res¬ 
ponse  of  the  specimens  was  measured  by  a  silicon  photodiode 
detector.**  These  two  detectors  were  placed  underneath  the  fur¬ 
nace  and  focused  through  a  sapphire  window  at  the  bottom  of  the 
furnace  onto  the  lower  surface  of  the  specimen.  It  should  be 
pointed  out  that  the  furnace  provided  the  ambient  temperature  of 
the  specimen.  The  mean  temperature  rise  of  the  specimen  due  to 
the  laser  flash  is  less  than  five  degrees. 

A  side-viewing  liquid  nitrogen  cooled  InSb  detector***  was 
used  for  measurements  of  thermal  diffusivity  at  room  temperature. 
For  this  purpose  the  specimens,  contained  in  a  thermally  insu¬ 
lated  holder,  were  placed  directly  in  the  path  of  the  laser  pulse 
and  immediately  in  front  of  the  detector  window.  For  measure¬ 
ments  at  all  temperatures,  the  specimens  were  coated  with  a  thin 


*  12-8722  InSb,  Barnes  Industry,  Stamford,  Conn. ,  USA. 

**  Built  for  this  purpose  at  VPI  &  SU  for  this  laboratory. 

***  Type  TI-ISV381  InSb,  Texas  Instruments,  Dallas,  Texas. 
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layer  of  graphite*  in  order  to  prevent  direct  transmission  of  the 
laser  flash  through  the  samples. 

The  signal  from  an  optical  detector  can  be  recorded  with  an 
appropriate  data  acquisition  system  such  as  the  storage  oscillo¬ 
scope**  used  for  the  present  study.  (A  computerized  data  acqui¬ 
sition  and  analysis  system  is  being  constructed  at  this  time.) 

The  oscilloscope  was  triggered  at  the  instant  the  laser  fired  by 
a  signal  from  a  photodiode  mounted  at  the  back  of  the  laser  from 
which  approximately  1%  of  the  total  flash  energy  is  emitted.  By 
appropriately  biasing  the  signal  from  the  detectors  corresponding 
to  the  temperature  of  the  specimen  prior  to  the  laser-flashing, 
only  the  transient  temperature  rise  is  recorded.  This  transient 
temperature  rise  is  shown  schematically  in  Figure  2.  The  perti¬ 
nent  value  to  be  measured  is  t1/2'  w^ich  represents  the  time 
required  for  the  lower  face  of  the  specimen  to  reach  one  half  of 
the  maximum  temperature  change.  The  absolute  temperature  rise 
need  not  be  known  quantitatively  for  the  determination  of  the 
thermal  diffusivity. 

For  unidirectional  heat  flow,  uniform  radiation  of  the  spe¬ 
cimen,  and  no  heat  losses  following  laser  flashing,  the  thermal 
diffusivity  (  k)  ,  the  half-rise  time  (t -^/2 )»  and  the  specimen 
thickness  (L)  are  related  by: 

k  =  1.37  L2A2t1/2  (1) 


*  dgf  123  Dry  Graphite  Film,  Miracle  Power  Products  Corp., 
Cleveland,  Ohio. 

♦♦Tektronix  5111  storage  oscilloscope. 
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At  higher  temperatures  where  heat  losses  may  be  appreciable 
during  the  time  of  the  specimen's  transient  temperature  response, 
equation  1  must  be  modified.  In  this  study  the  rate  of  specimen 


cooling  was  determined  by  measuring  the  relative  temperature 

rises  and  times  at  the  point  of  T _ ,  and  10  t,  Equation  1  then 

max  1/2 

was  readily  corrected  by  the  use  of  tables  compiled  by  Heckman 


(6).  The  precision  of  the  thermal  diffusivity  due  to  measurement 


errors  is  about  2.5%.  Data  obtained  on  Armco  iron  using  the 


equipment  of  this  study  was  accurate  to  within  4%  of  literature 
values  (7). 


The  thermal  conductivity  (X)  was  calculated  from: 


A  =  Kpcp  (2) 

where  p  and  Cp  are  the  density  and  specific  heat,  respectively. 

The  values  of  density  used  were  those  listed  in  Table  1,  whereas 
the  data  for  the  specific  heat  (1)  are  listed  in  Table  2.  Above 
room  temperature,  changes  in  the  specimen  thickness  and  density 
with  temperature  were  taken  into  account  in  the  evaluation  of  the 
thermal  diffusivity  and  conductivity. 

The  laser-flash  specimens  cut  from  a  given  refractory  compo¬ 
sition  were  measured  at  room  temperature.  A  total  of  four  sepa¬ 
rate  readings  of  t^^  were  taken  for  each  specimen  and  averaged. 
This  average  ^  ^/2  was  usec*  ^o  calculate  the  thermal  diffusivity 
of  the  specimen  at  room  temperature.  The  specimen  with  the  value 
of  thermal  diffusivity  closest  to  the  mean  value  of  all  the  spe¬ 
cimens  cut  from  the  same  refractory  was  judged  to  be  representa¬ 
tive  of  the  bulk  material  and  was  selected  for  measurement  of 
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thermal  diffusivity  as  a  function  of  temperature.  Equation  (2) 
was  then  used  to  calculate  the  thermal  conductivity  from  each 
thermal  diffusivity  data  point. 

2.2.2  Measurement  of  Thermal  Conductivity 

Values  for  the  thermal  conductivity  of  the  samples  provided 
by  manufacturers  A  and  B  were  obtained  by  the  guarded-plate  and 
hot-wire  methods,  respectively.  The  data  obtained  by  the  guarded 
plate  method  used  equipment  identical  to  that  designated  for  the 
ASTM  test  (4).  The  bricks  used  for  these  measurements  were  nomi¬ 
nally  identical  in  composition  and  microstructure  to  those  used 
for  the  laser-flash  measurements. 

Measurement  of  the  thermal  conductivity  by  the  hot-wire 
method  (3)  was  based  on  the  accurate  measurement  of  the  change  in 
resistance  of  a  pure  platinum  wire  connected  to  a  constant  cur¬ 
rent  source.  The  wire  is  located  between  two  22.5  cm  straights 
held  at  temperatures  ranging  from  room  temperature  to  1500° C.  A 
fully  computerized  data  acquisition  and  analysis  system*  was 
used.  Five  seperate  measurements  were  made  at  each  of  seven 
different  test  temperatures.  Specimens  for  the  laser-flash  mea¬ 
surements  were  cut  from  one  of  the  identical  bricks  used  for  the 
hot  wire  measurements. 


*  Hewlett-Packard  3052A  Data  Acquisition  System. 


12 


3.  RESULTS,  DISCUSSION  AND  CONCLUSIONS 

Tables  3  and  4  list  the  thermal  diffusivity  data  for  room 
temperature  measurements  of  the  individual  specimens  from  those 
refractories  which  exhibited  the  smallest  and  largest  specimen- 
to-specimen  variations,  respectively.  It  can  be  seen  from  the 
scatter  in  the  data  that  any  one  single  specimen  cannot  be  con¬ 
sidered  as  representative  of  the  bulk  material.  Accordingly,  in 
order  to  obtain  reliable  thermal  diffusivity  and  conductivity 
data,  care  must  be  taken  in  selecting  a  specimen  which  can  be 
considered  representative  of  the  bulk  material.  For  this  study, 
as  few  as  six  and  as  many  as  seventeen  individual  specimens  were 
cut  from  each  refractory  brick  for  each  direction  of  heat  flow 
where  appropriate.  Measurement  of  the  thermal  diffusivity  of 
twenty  specimens  at  room  temperature  requires  no  more  than  a  cou¬ 
ple  of  hours,  depending  on  the  number  of  measurements  made  for 
any  one  specimen.  Measurement  at  elevated  temperatures,  on  the 
other  hand,  is  more  time  consuming  and  thus  impractical  for  all 
specimens.  As  an  aside,  it  should  be  noted  that,  in  principle, 
specimen- to-specimen  variation  can  be  reduced  by  the  use  of  lar¬ 
ger  and  thicker  samples.  This,  however,  would  require  much  lar¬ 
ger  lasers,  with  prohibitive  cost  and  is  therefore  less  practi¬ 
cal  . 

Figures  3  to  13  show  the  experimental  data  for  the  thermal 
diffusivity  and  calculated  values  of  thermal  conductivity,  as  a 
function  of  temperature,  obtained  by  the  laser-flash  method  for 
all  eleven  refractory  compositions.  Also  included  in  these  fig- 
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ures  are  the  data  for  the  thermal  conductivity  obtained  by  the 
guarded-plate  or  hot-wire  method.  Comparison  of  the  data  shows 
that  at  elevated  temperatures,  with  the  exception  of  refractory 
4,  the  agreement  between  the  laser-flash  data  and  those  obtained 
by  the  guarded-plate  or  hot-wire  method  is  quite  good. 

For  refractory  4,  which  is  a  rather  porous  refractory,  the 
laser-flash  data  exceed  those  obtained  by  the  guarded-plate 
method  by  20  to  30%.  It  is  suggested  here  that  for  this  refrac¬ 
tory  the  laser-flash  data  are  too  high.  This  effect  most  likely 
arises  from  the  presence  of  the  pores  in  the  surface  of  the  spe¬ 
cimen.  These  pores  permit  the  absorption  of  the  laser  radiation 
within  the  interior  of  the  specimen  rather  than  at  a  planar  sur¬ 
face  assumed  in  the  development  of  the  underlying  theory.  The 
same  effect  arises  at  the  opposite  surface,  where  pores  permit 
detection  of  temperature  changes  in  regions  within  the  specimen. 
Both  effects  combine  to  reduce  the  effective  specimen  thickness 
which  tends  to  increase  the  value  of  the  thermal  diffusivity  cal¬ 
culated  from  t-jy2«  This  effect,  which  is  expected  to  be  a  func¬ 
tion  of  the  pore  volume  fraction  as  well  as  pore  size,  may  in 
principle,  be  accounted  for  theoretically.  However,  this  would 
require  information  on  the  surface  pore  configuration  for  any 
pertinent  material  to  be  tested,  coupled  with  a  transient  heat 
flow  analysis  of  the  spatially  non-uniform  absorption  of  the 
laser  radiation.  An  analysis  of  this  type  is  beyond  the  scope  of 
the  present  study. 


Near  room  temperature,  the  data  for  the  thermal  conductivity 
obtained  by  the  laser-flash  method  generally  exceed  those 
obtained  by  the  guarded-plate  and  hot-wire  methods.  It  is 
thought  that  this  effect  does  not  arise  from  some  systematic 
error  in  the  flash  method  because  good  agreement  was  found  bet¬ 
ween  data  obtained  by  this  method  and  data  for  thermal  conductiv¬ 
ity  standards  such  as  Armco  iron.  In  the  guarded-plate  method 
the  temperature  of  the  specimen  is  reported  in  terms  of  the 
linear  average  of  the  two  temperatures  recorded  by  the  thermocou¬ 
ples.  This,  however,  tends  to  over  and  under  estimate  the  spa- 
cial  average  of  the  temperature  and  thermal  conductivity,  respec¬ 
tively,  for  materials  with  a  negatively  temperature  dependent 
thermal  conductivity.  Possibly,  similar  effects  are  operative 
during  measurements  by  the  hot-wire  method. 

In  conclusion  the  results  of  this  study  indicate  that  for 
appropriately  chosen  specimens,  the  laser-flash  method  is  a  use¬ 
ful  tool  for  obtaining  the  thermal  diffusivity  and  conductivity 
of  dense  refractories. 
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TABLE  2.  SPECIFIC  HEAT  VALUES  USED  TO  DETERMINE 
THE  THERMAL  CONDUCTIVITY 


Specific  Heat  (cal/g.C) 


Temperature 


0 

93 

204 

316 

427 

538 

649 

760 

871 

982 

1093 

1204 

1316 

1427 

1538 

1649 

1760 


Refractory 

1,2 


0.193 

0.199 

0.206 

0.212 

0.220 

0.227 

0.234 

0.241 

0.248 

0.253 

0.258 

0.262 

0.266 

0.269 


Refractory 

3,10,11 


0.171 
0.196 
0.214 
0.226 
0.235 
0.242 
0.248 
0.252 
0.257 
0.260 
0.264 
0.267 
0.270 
0.  273 
0.276 
0.279 
0.282 


Refractory 

4 


0.184 

0.192 

0.214 

0.223 

0.229 

0.233 

0.237 

0.240 

0.242 

0.245 

0.247 

0.249 

0.251 

0.253 

0.255 

0.256 

0.258 


Refractory 

5-9 


0.208 

0.219 

0.232 

0.242 

0.251 

0.258 

0.263 

0.268 

0.273 

0.278 

0.283 

0.288 

0.293 

0.297 


Table  #  3 


Thermal  Diffusivity  of  Adjacent  Samples  in  Refractory  If  4 
at  Room  Temperature  Parallel  to  the  Pressing  Direction. 


Sample 

Thickness 

tl/2 

-  1 

Th.  Diff. 

_  3  ? 

Number 

(10  cm) 

(10  sec ) 

(10  cm“/sec) 

1 

2.11 

8.80 

7.19 

2 

2.33 

11.0 

6.85 

3 

2. 12 

8.50 

7.34 

4 

2.10 

8.15 

7.51 

5 

2.20 

9.50 

7.07 

6 

2.33 

10.5 

7.18 

7  * 

2.08 

8. AO 

7.15 

8 

2. 12 

8.80 

7.09 

AVERAGE  DIFFUSIVITY  -  0.00717  cm2 /sec 

AVERAGE  CONDUCTIVITY  -  0.764  W/m-K 

*  Denotes  specimen  used  for  high  temperature 
measurements 


Table  0  4 


Thermal  Diffusivity  of  Adjacent  Samples  in  Refractory  0  7 
at  Room  Temperature  Parallel  to  the  Pressing  Direction. 


Sample 

Number 

Thickness 

Cl0-1cm> 

»— « 

o 

1  ft 

M  ^ 

W  ^ 

(0  tJ 

o 

Th.  Di 

(10'2cm2 

1 

2.  11 

2.18 

2.83 

2 

2.23 

1.16 

5.95 

3 

2.22 

1.44 

4.75 

4 

2.20 

2.20 

3.05 

5* 

2.21 

- 

- 

6 

2.17 

1.38 

4.74 

7 

2.20 

2.00 

3.36 

8 

2.19 

1.75 

3.80 

9 

2.07 

1.30 

4.57 

10 

2. 12 

2.50 

2.50 

*  Specimen  damaged,  not  suitable  for  testing. 

(table  continued  on  next  page) 


Table  #4 
(continued) 

Thermal  Diffusivity  of  Adjacent  Samples  in  Refractory  #7 
at  Room  Temperature  Parallel  to  the  Pressing  Direction  . 


Sample 

Number 

Thickness 

(10  ^cm) 

Cl/2 

(10  ^sec ) 

Th.  Di 

(10  “cm^ 

11  ** 

2.18 

1.78 

3.71 

12 

2.22 

2.60 

2.63 

13 

2.11 

2.53 

2.44 

14 

2.23 

2.50 

2.76 

15 

2.26 

2.80 

2.53 

16 

2.22 

2.05 

3.3£ 

17* 

2.18 

- 

- 

18 

2.31 

1.95 

3.80 

19 

2.15 

1.50 

4.28 

AVERAGE 

AVERAGE 

DIFFUSIVITY  -  0.0359  cm2 

CONDUCTIVITY  -  9.60  W/m-K 

/sec 

*  Specimen  not  suitable  for  testing. 

**  Denotes  specimen  used  for  high  temperature 
measurements 
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Background  and  Objectives 
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In  several  previous  papers  [1-7]  the  present  authors  and  co-workers 
have  presented  the  results  of  studies  on  the  thermal  stresses  produced 
by  suddenly  applied  [1-5]  or  suddenly  removed  [6-7]  radiation  heating 
of  semi-transparent  material.  In  all  cases  we  assumed  well-collimated 
incident  radiation  striking  a  flat  plate  of  the  material  in  question 
and  being  absorbed  through-the-thickness  of  the  plate.  We  also  assumed 
normal  incidence  of  the  radiation,  and  ignored  internal  reflections.  In 
the  case  of  nuclear  radiation,  no  surface  reflections  occur  and  if  the 
material  is  either  purely  absorbing  or  optically  thin,  these  assump¬ 
tions  lead  to  a  realistic  model  of  the  heat  generation.  If  the  heating 
is  produced  by  thermal  radiation,  on  the  other  hand,  even  in  a  purely 
absorbing  medium  some  reflection  occurs  at  the  material  interfaces  due 
to  the  different  indices  of  refraction,  and  this  would  be  further  com¬ 
plicated  by  non-normal  incidence.  It  is  then  natural  to  ask  if  such 
reflections  significantly  alter  the  heat  generation  and  thermal  stress, 
and  what  the  effect  of  varying  the  angle  of  incidence  would  be.  In  the 
work  reported  here  we  provide  answers  to  some  of  these  questions. 

Our  objectives  here  are  to  determine  the  effects  of  interface  re¬ 
flections  on  the  thermal  stresses  generated  by  radiation  absorption  in 
a  semi-transparent  flat  plate,  and  to  study  the  effect  of  varying  the 
angle  of  incidence.  We  consider  both  suddenly  applied  radiation,  and 
a  sudden  interruption  of  radiation  from  an  established  steady-state 
condition.  Reflections  at  the  surface  of  the  plate  should  result  in 
"trapping"  more  of  the  radiation  within  the  plate,  and  thus  to  a  larger 
volumetric  heat  source.  Whether  or  not  this  leads  to  increased  thermal 
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stress  can  only  be  determined  by  analysis.  However,  varying  the  angle 
of  incidence  away  from  normal  increases  the  surface  reflection  from 
the  plate,  but  it  also  increases  the  path  length  of  the  radiation  with¬ 
in  the  plate,  changes  which  should  have  opposite  effects  on  the  heat 
generation.  Again,  only  analysis  can  determine  the  net  result,  as  we 
shall  see. 


Method  of  Analysis 

Consider  a  flat  plate  of  thickness  w  =  2a  (See  Fig.  1)  insulated 
on  the  front  face  (x  =  0)  and  cooled  on  the  rear  face  (x  =  w)  by  a  con¬ 
vective  medium  at  temperature  T  and  heat  transfer  coefficient  h.  The 
plate  is  assumed  to  be  initially  in  equilibrium  at  temperature  Tot;  at 
time  t  *  U,  a  collimated  beam  of  thermal  radiation  of  intensity  I  is 
applied  to  the  surface  at  x  ■  0  at  an  incidence  angle  0.  This  radiation 

originates  in  a  medium  of  refractive  index  n  ,  while  the  plate  has 

o 

index  of  refraction  n^  and  the  convective  medium  .  To  determine  the 

thermal  stresses  generated  by  absorption  of  the  radiation,  we  must  (a) 

formulate  the  appropriate  volumetric-heat-source  term  for  the  conduction 

equation;  (b)  solve  the  conduction  equation  for  the  unsteady  temperature 

profile  T(x,t);  and  (c)  compute  the  thermal  stress  o  (x.t)  from  the 

yz 

temperature  profile.  Since  we  want  to  compute  the  thermal  stresses  both 
for  suddenly  applied  and  suduenly  removed  incident  radiation,  the  analysis 
in  parts  (b)  and  (c)  must  be  performed  for  both  of  these  cases. 

We  first  derive  an  expression  for  the  heat  source  g'"(x,<f>)  by  summing 
the  contributions  from  each  pass  of  the  radiation  through  the  nlate 
resulting  from  successive  reflections  at  the  plate  surfaces,  as  shown 
in  Fig.  2.  We  assume  a  purely  absorbing  medium  with  absorption  (ex¬ 
tinction)  coefficient  y.  The  angle  of  incidence  <f  and  the  refracted 
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angle  8  are  related  through  the  complex  refractive  indices  of  the  two 
materials  [8].  At  a  given  depth  x  into  the  medium,  it  is  necessary  to 
add  heat  generation  due  to  radiation  absorption  at  points  1,  2,  3,  ..  . 
A  fraction  is  reflected  at  the  first  interface,  so  that  at  point  1, 
the  radiation  intensity  is 

I0  ^  “  pl^  exp  M  Ll^ ’ 

where 


=  x/cosfi  .  (1) 

Similarly  at  point  2,  the  intensity  is 

IQ  ^  “  Pl^  p2  £Xp  +  L2^’ 

where 

*  (w  -  x)/cosO,  (2) 

and 


L  «  +  L  . 


Thus,  in  general,  the  intensity  at  x  is 


(3) 


Summing  the  geometric  series,  we  find 


I 


a 


4 


I 


1 


I(x)  =  I0  (1  -  I'.p 


r  ",jL1 

e  +  p.^e 


-P(L  +  L0) 


-2pL 


1  -  PlP2e 


(5) 


or  explicitly  in  terms  of  x,  w,  and  6, 


I(x)  =  IQ(1  -  px) 


r  -px/cos0  -p(2w  -  x)/cos0 

e  +  o2e 


-2pw/ cosO 


1  -  P1P2e 


(6) 


The  volumetric  heat  source  at  location  x  is  then 


g'"(x,.t.)  =  pl(x). 


(7) 


i.e.  the  rate  of  absorption  of  the  radiation  at  x  with  I(x)  given  by  Eq, 
(6).  It  must  be  recognized  in  Eq.  (6)  that  the  reflectivities  p^  and 
f * 2  also  depend  on  (and  thus  4>)  through  the  refractive  indexes  n^,  n^ , 
and  n2  [ 8 ] . 


(a)  Suddenly  applied  radiation 


We  seek  a  solution  of  the  heat  conduction  equation 

k  (x,t)  +  g'"(x)  =  pc  ,  t  >  0, 

0x 


(8) 


where  T(x,t)  represents  the  temperature  at  location  x  at  time  t,  k  the 
thermal  conductivity,  p  the  density,  and  c  the  specific  heat  capacity. 
The  appropriate  boundary  conditions  for  this  case  are 


rr 

3x 


(0,1)  =  0  , 


(9) 
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(10) 


-  Ic|“  (w,t)  =  h  [T(w,  t)  -  TJ, 


with  initial  condition 


T (x , 0)  =  T^. 

Using  standard  techniques  [9]  we  find  the  solution  to  be 


(11) 


T1 (x,t)  =  T  +  l  B  (t)  cos  (X  x),  0  <  x  <  w,  t  >  0,  (12) 

i  00  n  n  n  —  — 

n=0 


where 


Bn(t) 


-<X^t  2 

(filn  +  *2n>  (1  -  e  >/Xn  V 


(13) 


In  Eq.  (13), 


sin  (2A  w) 

N  =  —  +  - S — 

n  2  4A 

n 


(14) 


n 


In  2  2  2 . 

p  +  A  cos  6 
n 


-pw/ COS0 

1 _ |e  [-pcos0cos(A^w)  +  X_  sin(A_w) 


n  n  n 


(15) 


2 

COS  0  ]  +  p  COS  0 |, 


u 


'2 


pw/ COS0 


2n  2  ,  ,2  2 
p  +  A  cos  0 
n 


[pcos0cos(X  w)  +  X  sin  (X  w) 
n  n  n 


2 

COS  0  ]  -  pCOS  <*  |  , 


(16) 
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where 


1 

k 


1 


■"o  (1  -  ^ _ 

P ^ > 2  exP  (~2jiw'/cos0) 


(17) 


and 


<\>2  ~  >'o  exp  (-2pw/cosf)  .  (18) 

The  eigenvalues  A^  are  the  roots  of  the  transcendental  equation 

A  tan  (A  w)  =  —  .  (19) 

n  n  k 

The  thermal  stresses  are  computed  from  the  well-known  equation  for  an 
unconstrained  flat  plate  [10], 


V  f-T(x-t)  +  8  +  fr 

J  L  -a  2a 

.  ;  cT(£;,t)df>] 

-a  J 


(20) 


where,  for  convenience,  we  have  now  shifted  the  origin  of  coordinates  to 
the  center  of  the  plate  such  that  -  a  <_  x  a.  Substituting  Eq.  (12) 
into  Eq .  (20)  and  performing  the  indicated  integrations,  we  find 


a  (x,t) 
yz 


«E_ 

1-v 


y  b 

n 


(t) 


F 

n 


(x)  , 


(21) 


where 


F  (x)  =  -cos  X  (x  +  a)  4 
n  n 


sin  2X  a  ,  /  2  sin  X  a 

n  3x  /  n 


2a3  X2 
\  n 


a  sin  2X  a 
n 


,  -a<x<a,  t>0. 


We  will  study  this  result  in  more  detail  in  a  future  section. 


(b)  Suddenly  interrupted  radiation 


For  generality,  we  assume  that  the  incident  radiation  is  interrupted 
after  a  time  interval  x ;  x  -*■  °°  corresponds  to  equilibrium  conditions 
and  can  easily  be  recovered  from  the  following  general  result.  We  thus 
seek  a  solution  to  Eq.  (8)  subject  to  boundary  conditions  (9)  and  (10), 
but  with  initial  condition 


T(x,0)  =  4  l  (x)  cos  Xn  x,  0  x  <_  w  .  ^3) 


The  solution  is  easily  shown  to  be 


,  2 

00  -k  X  t 

T(x,t)  =  T  4  I  B  (t)  e  n  cos  X  x,  0  <  x  <  w  , 
<x.  n  n  —  — 

n=0 


with  corresponding  thermal  stress  distribution 


->:X2t 

°yz  -1=7  EnBn  6  "  Fn  (x)  ’ 

n=U 


-  a  <  x  <  a 
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In  the  following  section  we  display  computed  results  from  these 
equations,  and  in  particular,  we  examine  the  effect  of  varying  the 
incidence  angle  <j>  and  the  role  of  multiple  surface  reflections.  It  is 
clear  from  Eqs.  (15-19)  that  4'  enters  in  a  rather  complicated  fashion, 
particularly  since  and  p ^  also  depend  on  0  and  4*  through  the  complex 
refractive  indices  for  absorbing  media  [8], 
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Numerical  Results 


To  illustrate  the  effects  of  interest  to  this  investigation  we  use 
a  particular  example  of  a  glass  '-’ith  wavelength  dependence  of  trans¬ 
mittance  as  shown  in  Fig.  (3)  [11].  This  introduces  another  variable 
to  be  considered,  but  it  does  help  show  the  role  of  the  extinction 
coefficient  m,  since  we  exhibit  results  fur  2.6  and  3.5  pm,  respectively, 
where  the  absorption  coefficients  are  quite  different. 

In  Figs.  (4)  and  (5)  we  display  the  non-dimensional  temperature 
and  the  stress  distributions  as  computed  from  Eqs.  (12)  and  (22)  im¬ 
mediately  after  the  radiation  begins.  The  non-dimensional  temperature 
is  defined  as 

T(x,t)  -  Tm 

T*(x,t)  =  --T— Tf -  (26) 

F  00 

where  T  is  the  temperature  at  which  the  material  fails.  It  is  clear 
F 

that  the  greater  temperature  rise  occurs  at  normal  incidence  (<J>  =  0°). 
However,  the  stresses  are  greatest  for  a  60°  incidence  angle.  This  is 
illustrated  more  clearly  in  Fig.  (6),  where  we  plot  the  maximum  tensile 
stress  as  a  function  of  incidence  angle  for  two  different  wavelengths. 
Both  curves  show  a  broad  maximum  in  the  region  near  <j>  =  60°.  Thus, 
apparently  at  this  incidence  angle,  the  competing  effects  of  increased 
reflection  at  the  surface  and  greater  path  length  in  the  plate  are  just 
balanced.  The  overall  stress  levels  are  higher  at  2.b  ym  because  the 
absorption  coefficient  p  is  much  larger  at  this  wavelength. 

We  would  like  to  assess  the  effect  of  ignoring  multiple  surface 
reflections  on  the  temperature  and  stress  profiles.  In  Fig.  (7)  we 
have  plotted  the  non-dimensional  temperature  profile  immediately  after 
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application  of  the  radiation  at  normal  incidence.  The  two  curves  re¬ 
present  temperature  profiles  resulting  from  the  heat  source  distribution 
g"'(x,0°),  which  includes  reflections,  and  the  heat  source  g"'(x),  which 
does  not.  Figure  (8)  shows  the  corresponding  stress  distributions.  It 
is  seen  that  the  temperature  profile  shows  a  miximum  difference  of  about 
20%  at  the  right  face,  away  from  the  incident  radiation.  The  maximum 
difference  in  the  stress  distributions  occurs  at  the  center  of  the  plate, 
and  is  of  the  order  of  5%. 

The  stress  levels  predicted  in  our  analysis  are  still  well  below 

the  failure  stress  for  the  material  considered  [11]  because  this  material 

has  a  very  low  thermal  expansion  coefficient.  A  radiation  source  of  the 

6  2 

intensity  used  in  these  calculations  (5  x  10  W/cm  )  would  melt  the  mate¬ 
rial  before  it  would  fracture.  Many  other  brittle  materials,  however, 
have  much  larger  values  of  thermal  expansion  coefficient  and  the  stresses 
will  be  correspondingly  higher.  In  these  cases,  the  structure  could 
will  fail  by  thermal  stress  fracture. 

If  the  radiation  is  incident  for  a  time  interval  of  sufficient 
length  to  allow  an  equilibrium  temperature  distribution  to  be  established, 
and  the  heat  source  is  then  suddenly  removed,  the  rapid  collapse  of  the 
temperature  profile  can  again  produce  large  stresses.  The  resulting 
stress  distribution  will  depend  on  the  angle  of  incidence  of  the  radia¬ 
tion  through  the  initial  temperature  profile.  This  is  illustrated  in 
Fig.  (9)  where  we  have  plotted  the  maximum  tensile  stress  at  several 
times  during  the  cooling  phase  vs.  angle  of  incidence.  It  is  seen  that 
at  early  times  the  maximum  stresses  still  occur  for  an  incidence  angle 
near  60°,  but  as  time  passes  this  effect  is  gradually  lost.  After 
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Dut  10  seconds,  for  this  case,  all  history  of  the  incidence  angle  has 

sappeared . 

The  evolution  of  the  temperature  and  stress  profiles  during  the  cool- 
m  phase  is  clearly  apparent  from  Figs.  (10)  and  (11).  In  these  figures 
show  temperature  and  stress  distributions  at  several  times  after  the 
at  source  is  removed,  with  time  displayed  on  a  third  axis.  With  the 
ft  face  no  longer  receiving  radiation,  the  temperature  there  drops  most 
aidly,  and  the  profile  changes  from  concave-up  to  concave-down.  The 
ress  distribution  (Fig.  11)  thus  shifts  from  concave-down  to  concave- 
.  At  intermediate  times,  however,  the  stress  distribution  shows  two 
:al  minima  which  can  be  correlated  with  inflection  points  in  the 
aperature  profiles.  Thus  it  appeals  that  the  stress  profile  always 
splays  concavity  opposite  to  that  of  the  temperature  profile. 

lclusions 

It  is  clear  from  the  above  that  multiple  surface  reflections  and 
^-normal  incidence  angles  can  be  accounted  for  quite  straight- 
•ewardly  in  studies  of  thermal  stresses  generated  by  radiation  heat- 
;.  Temperatures  throughout  the  medium  decrease  monotonically  as 
angle  oi  incidence  increases,  whereas  the  largest  thermal  stress 
ses  through  a  broad  maximum  near  an  incidence  angle  of  bU°.  Thus 
would  be  advantageous  in  design  of  structures  exposed  to  transient 
iation  heating  to  avoid  incidence  angles  near  60°.  Upon  sudden  inter- 
tion  of  the  radiation,  this  same  effect  appears  during  the  early 
ges  of  the  cooldown,  but  gradually  disappears  as  the  history  of 
original  heating  pattern  ir  lost. 


For  radiation  at  normal  incidence,  ignoring  surface  reflections 


leads  to  a  maximum  error  of  approximately  +  20%  in  the  temperature  dis¬ 
tribution,  but  no  more  than  +  10%  in  the  thermal  stress. 
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Figure  2.  Ray  tracing  to  determine  the  total 
energy  absorption  rate 


Figure 
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3.  Wavelength  dependence  ot  transmittance  for 
1  cm.  thick  Fused  Silica  I! 7940  glass  plate 
20°C  (from  Ref .  11) • 
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Figure  4.  Non-dimensional  temperature  distribution  for 

the  heated  plate  at  t  =  10~6  seconds  for  various 
incidence  angles  <j>. 


Figure  5,  Stress  distribution  for  the  heated  plate  at  10 
seconds  for  various  incidence  angles. 
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Figure  8.  Effect  of  interface  reflections  on  the  stress 
distribution  in  a  heated  flat  plate  at  t  =  10 
seconds,  normal  incidence,  A  =  2.6  pm. 
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Figure  9.  Maximum  tensile  stress  at  various  times  during 
the  cooling  of  the  radiatively  heated  plate  as 
a  function  of  the  incidence  angle  of  the  radia' 
Lion,  A  =  2.6  pm. 
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Figure  10.  Three-dimensional  portrayal  of  the  non-dimensional 
temperature  distribution  during  the  cooling  of  the 
radiatively  heated  plate,  normal  incidence,  A  =  2.6  um. 
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Figure  11.  Three-dimensional  portrayal  of  the  sir  ss  distribution 
during  the  cooling  of  the  radiatively  heated  plate, 
normal  incidence,  A  =  2.6  pm. 
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ROLE  OF  ABSORPTION  COEFFICIENT  IN  THE  FREQUENCY  DEPENDENCE  OF  THE  THERMAL 
STRESSES  IN  A  PARTIALLY- ABSORBING  PLATE  SUBJECTED  TO  CYCLIC  THERMAL  RADIATION 
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ABSTRACT 

The  frequency  dependence  of  the  thermal  stresses  in  a  partially  absorb¬ 
ing  plate  subjected  to  cyclic  thermal  radiation  is  shown  to  be  a  function 
of  the  optical  thickness. 


INTRODUCTION 

Earlier  reports  in  this  series  of  analytical  studies  of  the  thermal 
stresses  in  partially  absorbing  plates  subjected  to  radiation  heating 
concentrated  on  a  series  of  boundary  conditions  including  symmetric  [1,2] 
or  asymmetric  steady-state  radiation  [3,4]  combined  with  convective 
cooling,  sudden  loss  of  radiation  following  long-term  steady-state 
radiation  [5,6]  and  cyclic  radiation  combined  with  convective  cooling  [8,9]. 
These  latter  two  studies  concentrated  on  the  frequency  dependence  of  the 
thermal  stresses  for  a  single  value  of  optical  thickness  near  the 
value  which  under  steady-state  radiation  resulted  in  the  highest  value 
of  thermal  stress. 

The  purpose  of  this  paper  is  to  present  an  analysis  and  numerical 
result  for  the  effect  of  optical  thickness  on  the  frequency  dependence  of 
the  thermal  stresses  in  a  partially  absorbing  plate  subjected  to  cyclic 
radiation. 

ANALYSIS 

Boundary  and  Initial  Conditions 

The  flat  plate,  infinite  in  extent,  is  located  in  the  yz-plane  and 
-a<x<a.  At  time  t-0  the  front  face  of  the  plate  (x  *  -a)  is  subjected 
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to  normally  incident  spatially  uniform  cyclic  radiation  and  cooled  by 
Newtonian  convection.  The  back  face  of  the  plate  (x  *  a)  is  fully 
insulated.  The  reflectivity  of  the  material  of  the  plate  is  assumed 
to  be  sufficiently  low  that  multiple  reflections  within  the  plate  can 
be  neglected.  For  the  effect  of  non-normal  incidence  and  multiple 
reflections  under  conditions  of  steady-state  radiation  the  reader  is 
referred  to  the  study  of  Thomas  et  al.  [7].  The  material  properties 
which  affect  the  temperature  and  stress  such  as  the  emissivity  (absorp¬ 
tivity)  ,  absorption  coefficient,  coefficient  of  thermal  expansion. 

Young's  modulus  of  elasticity,  Poisson's  ratio,  thermal  conductivity 
and  the  thermal  diffusivity  are  assumed  to  be  independent  of  temperature. 
The  optical  properties  of  the  material  of  the  plate  were  assumed  to  be 
"grey",  i.e.  independent  of  wavelength.  The  temperature  of  the  plate 
over  its  total  thermal  history  is  assumed  to  be  sufficiently  low  that 
the  intensity  of  the  re-emitted  radiation  can  be  considered  negligible 
compared  to  the  magnitude  of  the  incident  radiative  heat  flux. 

The  flux  q(t)  of  the  normally  incident  radiation  at  it  *  -a 
was  taken  to  be  of  the  form: 

q(t)  -  qQ[l  +  6  sin(wt)]  (1) 

where  qQ  is  the  value  of  heat  flux  at  t  *  0,  u  is  the  angular  frequency 

and  0  <  6  <  1.  For  the  numerical  example  to  be  presented  later  6  was  taken 

to  be  identically  equal  to  unity,  so  that  the  intenrity  of  the  incident 

radiation  varies  between  0  and  2q  . 

o 

The  heat  flux  q  at  any  position  x  within  the  plate  due  to  thermal 
radiation  entering  at  x  =  -a  is: 

q(x,t)  *  eqQe  U^a+X^ [i  +  5  sin  (wt) ]  (2) 

where  e  ■  1  -  r  with  r  being  the  reflectivity  and  y  the  absorption 
coefficient.  The  rate  of  internal  heat  absorption  per  unit  volume  (g'"  ) 
within  the  plate  at  x  is : 

g"' (x,t)  -  yeqQe“y^a+x^  [1  +  <5  sin(ojt)  ]  (3) 

Solutions  for  the  transient  temperature  T(x,t)  are  obtained  by  solving 
the  heat  conduction  equation  [10]: 


2 


(4) 


n  g1"  (x.t)  a  1  3T 
3x2  k  K  3t 

where  k  is  the  thermal  conductivity  and  k  is  the  thermal  diffusivity. 

For  the  given  heating  and  cooling  conditions,  the  initial  and 
boundary  conditions  are: 


T(x,0)  -  T 


3T  , 

3x  (a,t) 


(5) 

(6) 

(7) 


Solutions  for  the  transient  temperatures  were  obtained  using  standard 
techniques  [ 10 ] . 


T(x, t)  *  T  +  l  B  (t)  cos (A  x) 
o  .  n  n 

n-1 

where  the  A^  are  the  roots  of  the  transcendental  equation 


(8) 


A  tan(2A  a)  *  — 
n  n  k 


(9) 


and 


G  5 


B  (t) 


»  +  (<A/) 


2  -icA  t 

_  9  —  [kA  sin(o)t)  -  u  cos  wt+ioe  n  ] 

4  n 


G  ,2„ 

n  ,,  -icA  t, 
+  — 5 —  (1-e  n  ) 

<A 

n 


(10) 


with 


kN 


A  sin(2A  a)+y  cos(2A  a)-ye 
n  n  n 


-2ya 


2  2 
y  +  A 

n 


(ID 


and 


N  •  a  +  — cos2(2A  a) 

n  2kA2  n 

n 


(12) 


THERMAL  STRESSES 


Thermal  stresses  a  (x,t)  are  obtained  by  substituting  for  T  from 
y ,  z 

equation  (8)  into  [11] 
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the  non-dimensional  stress 

*  ,  cr(l-v)k 
aEqQ£  a 

*  ,  2 

the  non-dimensional  time,  t  =  ict/a  and  the  non-dimensional  distance, 
x/a. 

NUMERICAL  RESULTS  AND  DISCUSSION 

The  analytical  results  for  the  temperatures  and  stresses  obtained 
above  can  be  demonstrated  most  conveniently  by  means  of  numerical 
examples.  With  the  exception  of  numerical  values  for  the  optical 
thickness,  ya,  and  the  angular  frequency,  w,  for  convenience  the  other 
relevant  parameters  were  chosen  to  have  constant  values  as  follows: 

heat  transfer  coefficient  h  *=  0.006  W.cm  °C  ;  thermal  conductivity 

-1  -1  2  -1 
k  *  0.3  W  cm  °C  ;  thermal  diffusivity  <  =  0.1  cm  s  ;  and  plate 

half-thickness  a  =  1  cm.  The  amplitude  factor,  6,  is  assumed  equal  to  be 

unity  throughout. 
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For  the  specific  case  of  pa  =  3  and  <*>  =  0.21  rad.s  \  Figs.  1  and 
2  show  the  temperatures  and  stresses,  respectively,  at  the  front 
face  (x  =  -  1  cm) ,  center  (x  =  0  cm)  and  back  face  (x  =  1  cm)  for  a 
number  of  cycles  corresponding  to:  (a)  the  initiation  of  radiation  at 
t  »  0;  (b)  after  a  number  of  cycles  at  which  the  plate  has  reached  an 
intermediate  level  of  temperature  and,  (c)  after  sufficient  cycles  that 
no  signicant  cycle-to-cycle  variation  in  the  magnitude  of  temperature 
and  thermal  stresses  occurs.  For  all  three  cases,  the  cyclic  variation  of 
the  temperature  decreases  from  the  front  to  the  back  face.  Of  course,  this 
is  expected  from  the  specific  boundary  conditions  chosen  for  this  analysis. 
Also,  as  expected,  the  cyclic  temperature  variation  at  the  center  and  back 
face  lag  in  time  behind  that  at  the  front  face.  Also  note  that  for  this 
value  of  angular  frequency  (0.21  rad.s  ^)  for  the  data  in  Figs.  1  and  2,  the 
cyclic  variation  of  temperature  within  any  one  cycle  is  small  compared  to  the 
mean  temperature  rise  of  the  plate. 

The  magnitude  of  thermal  stresses  at  the  front  and  back  face  become 
increasingly  more  positive  (i.e.,  tensile)  as  the  plate  is  subjected  to  an 
increasing  number  of  cycles,  whereas  the  stresses  at  the  center  become 
increasingly  more  negative.  This  latter  effect  arises  because  as  time 
progresses  the  plate  develops  a  temperature  distribution  which  is 
generally  concave  downward  on  which  the  cyclic  variation  is  superimposed. 

As  pointed  out  elsewhere  [12]  such  a  concave  downward  temperature  distribu¬ 
tion  will  result  in  thermal  stresses  which  are  tensile  in  the  front  and 
back  surface  of  the  plate,  in  spite  of  the  fact  that  after  long  periods  of 
time  the  mean  temperatures  at  the  back  surface  and  center  are  higher  than 
those  in  front. 

The  magnitude  of  the  thermal  stresses  shown  in  Fig.  2  can  be  compa'/ed 
with  those  obtained  in  the  earlier  study  for  identical  values  of  all  relevant 
parameters,  but  with  convective  cooling  of  the  plate  at  the  back  face  [8]. 

It  is  easily  noted  that  the  values  of  the  stresses  for  front  or  back  fare 
cooling  are  almost  identical.  Such  near  equality  for  the  stresses  for  front 
and  back  face  cooling  was  also  noted  for  steady-state  radiation  [3,4]. 

More  recently  it  was  pointed  out  [5,6]  that  as  the  plate  thickness  a  ■*  0 , 
the  stresses  for  convective  cooling  at  the  front  or  back  become  identical. 
For  a  ~  10  cm,  differences  were  noted  for  these  cooling  conditions  for 
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in  a  plate  due  to  the  sudden  interruption  of  long-term  steady- 
diation  [5,6]. 

-4  -1 

a  very  low  value  of  angular  frequency,  w  =  10  rad.s  ,  and  pa  *  3, 

and  3b  show  the  temperature  and  thermal  stresses,  respectively, 

ction  of  time  over  the  duration  of  the  first  cycle  from  the 

of  the  initiation  of  radiation  at  t  =  0  at  a  value  q  =  qQ  to  the 

of  time  within  the  cycle  at  which  the  magnitude  of  heat  flux  q  =  0. 

3a  the  magnitude  of  the  temperatures  at  the  front,  center  and  back 

close  together  to  be  plotted  separately  for  the  relatively  large 

ed  for  the  temperature  axis.  The  time  at  which  the  tempe~ature 

its  maximum  corresponds  exactly  to  the  time  at  which  the  heat  flux 

a  value  of  2qQ.  This  Implies  that  at  sufficiently  low  frequenc  js, 

ed  out  also  in  an  earlier  study  [9],  that  following  a  brief  initial 

t,  the  temperature  and  stresses  in  the  plate  cycle  along  with  the 

neous  value  of  heat  flux,  with  a  corresponding  value  for  a  constant 

k  equal  to  the  instantaneous  value  of  cyclic  heat  flux.  Confirmation 

conclusion  is  given  by  Fig.  3b  which  shows  that  after  a  brief 

2  2  3 

transient  which  is  complete  at  <t/a  *  2x10  (or  2x10  s) ,  the  stresses 

:tly  proportional  i__  the  instantaneous  heat  flux.  Furthermore,  the 

ie  of  stress  corresponds  to  twice  the  value  found  for  the  steady- 

:ess  encountered  for  steady-state  radiation  of  magnitude  qQ  [3,4]. 

ire  4  shows  the  magnitude  of  the  thermal  stress  in  the  front  face 

-4  -1 

Identical  value  of  frequency  (10  rad.s  )  and  time  period  within 
:  cycle  as  for  the  data  shown  in  Fig.  3b,  but  for  a  range  of  values 
il  thickness.  The  value  of  peak  stress  for  any  value  of  pa  is 
*  value  for  radiation  of  constant  intensity  [3,4];  again,  this 
the  earlier  statement  that  at  low  frequency,  for  purposes  of  evaluating 
lerotures  and  stresses,  the  plate  can  be  regarded  to  be  subjected  to 
state  heat  flux. 

re  5  shows  the  frequency  dependence  of  the  maximum  stresses  during 
t  the  front  face  (x  *  -a)  after  a  sufficient  time  period  such 
ycle-to-cycle  variation  in  the  magnitude  of  temperature  or  stresses 
These  results  clearly  indicate  that  the  relative  frequency  dependent 
ermal  stresses  is  a  function  of  the  optical  thickness.  As  the 
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frequency  u>  -*■  »  the  magnitude  of  stress  approaches  the  value  for 
the  stresses  sfter  long-term  radiation  with  intensity  qQ.  As  u»  0 
(but  not  u>  •  0)  the  magnitude  of  maximum  stress  within  the  cycle 
approaches  the  value  for  long-term  steady-state  radiation  of  intensity 
2qQ.  At  to  =  0,  the  intensity  of  radiation  remains  constant  at  a  value 
qQ  with  stresses  corresponding  to  steady-state  radiation  of  intensity 
qQ.  For  this  reason,  the  magnitude  of  stress  displays  a  discontinuity 
at  <ii  *  0,  with  a  value  of  stress  as  u>  -*■  0,  which  is  twice  the  value 
for  oj  »  0. 
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NON-DIMENSIONAL  TIME,  <1/0* 


~ dimensional  time,  *t/o* 


NON-HMCMHONAl  TIMC,  **/«* 

Transient  temperature  (a)  and  thermal  stresses  (b)  in 
partially  absorbing  flat  plate  heated  in  front  by 
normally  incident  cyclic  thermal  radiation  and  cooled 
by  convection  for  a  value  of  angular  frequency  u  ■ 
ICTUS'*!  and  values  of  pertinent  parameters  identical 
to  those  of  Figs.  1  and  2.  (Note,  temperatures  in 
the  various  positions  in  the  plate  vere  too  close 
to  be  plotted  separately  for  the  scale  used  for  Fig. 
3a.) 


NON*  DIMENSIONAL  TIME.  <f/«* 

Fig.  4.  Maximum  thermal  stresses  In  front  surface  of  partially  absorbing 
flat  plate  heated  In  frost  by  normally  incident  cyclic  thermal 
radiation  and  cooled  by  convection  for  an  angular  frequency 
u  *  and  range  of  values  of  optical  thickness  and  other 

parameters  identical  to  those  for  Figs.  1,  2  and  3. 


Fig.  3.  Frequency  dependence  of  maximum  thermal  stresses  in  front  surface  of 
partially  absorbing  flat  plate  heated  in  front  by  normally  Incident 
cyclic  radiation  and  cooled  By  convection  for  range  of  optical  thickness 
and  other  parameters  identical  to  those  of  previous  figures. 
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CHAPTER  VI 


THERMAL  STRESSES  IN  PARTIALLY  ABSORBING  FLAT  PLATE  DUE 
TO  SUDDEN  INTERRUPTION  OF  STEADY-STATE  ASYMMETRIC  RADIATION,  I: 
CONVECTIVE  COOLING  AT  REAR  SURFACE 


by 


I.  Frankel,  J.  P.  Singh,  J.  R.  Thomas,  Jr.,  D.  P.  H.  Hasselman 


Departments  of  Mechanical  and  Materials  Engineering 
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ABSTRACT 


An  analysis  is  presented  for  the  thermal  stresses  in  a  partially 
absorbing  flat  plate  due  to  sudden  interruption  of  steady-state  asymmetric 
thermal  radiation  on  the  front  face  with  convection  cooling  at  the  rear  face. 
Maximum  transient  tensile  stresses  occur  in  the  rear  surface  and  are  almost 
independent  of  the  heat  transfer  coefficient,  h.  For  optical  thickness 
pa  >  0.5,  the  maximum  tensile  stress  initially  rises  to  reach  a  peak  and 
then  decreases  with  time.  In  contrast,  for  pa  <  0.5,  the  stresses  decrease 
monotonically.  It  is  shown  that  a  plot  of  non-dimensional  stress  as  a  func¬ 
tion  of  non-dimensional  time  for  a  given  value  of  pa  and  h  does  not  result 
in  a  universal  plot  from  which  thermal  stresses  can  be  obtained  for  any  value 


of  plate  thickness. 


INTRODUCTION 


Solutions  for  the  transient  temperatures  and  thermal  stresses  for 
heating  frequently  can  be  obtained  from  the  corresponding  solutions 
for  cooling  (or  vice-versa)  by  a  simple  transformation  and/or  change 
in  sign.  This  generally  will  be  the  case  if  the  mode  of  heat  transfer, 
such  as  convection,  is  the  same  for  heating  and  cooling. 

In  practice,  however,  conditions  are  encountered  for  which  the 
heat  transfer  conditions  for  heating  and  cooling  are  different.  This 
will  require  separate  solutions  for  the  transient  temperatures  and 
corresponding  thermal  stresses.  Such  solutions  are  expected  to  be 
coupled,  as  the  initial  and  boundary  conditions  for  cooling  following 
heating  (or  vice-versa)  are  governed  by  the  magnitude  and  distribution 
of  the  temperatures  and  thermal  stresses  reached  during  heating. 

The  above  situation  is  encountered  in  such  installations  as 
concentrated  solar  energy  receivers  or  any  other  structures  or  components 
subject  to  simultaneous  radiation  heating  and  convective  cooling.  Such 
a  structure  may  have  reached  a  condition  of  steady-state  temperature  and 
thermal  stress,  when  the  radiation  is  suddenly  decreased  to  zero,  leaving 
the  structure  or  component  subject  to  rapid  cooling  by  convection  only. 

Such  an  event  is  expected  to  be  a  common  occurrence  for  solar  energy 
collectors  under  conditions  of  intermittent  cloud  cover. 

The  purpose  of  this  paper  is  to  illustrate  such  coupling  between 
the  thermal  stresses  encountered  during  heating  and  cooling.  Specifically, 
thermal  stresses  will  be  analyzed  for  a  semi-absorbir  •»  flat  plate  initially 
at  steady-state  conditions  after  being  subjw«.  c  d  ;  a  long-term  constant 
asymmetric  radiation  at  the  front  face  for  which  the  temperature  and 


stresses  were  obtained  in  earlier  studies  [l,2].  The  radiation  is 
suddenly  removed,  subjecting  the  plate  to  cooling  by  convection  at  either 
the  rear  or  front  face.  For  convenience  of  the  reader.  Part  I  will  con¬ 
centrate  on  cooling  at  the  rear  face,  with  the  solutions  for  front  face 
cooling  presented  in  Part  II. 

ANALYSIS 

Initial  and  Boundary  Conditions 

The  geometry  selected  for  the  analysis  consisted  of  an  infinite  flat 
plate  located  in  the  y-z  plane,  with  -a<x<a.  The  plate,  initially  at  a 
uniform  temperature  Tq,  is  subjected  to  a  normally  incident  uniform  radia¬ 
tion  qQ  and  cooled  by  convection  at  the  rear  surface  (x  =  a) .  When  the 
steady  state  conditions  for  temperatures  and  stresses  inside  the  plate  are 
reached,  the  radiation  flux  qQ  is  discontinued.  The  temperature  solution 
for  the  above  st  , Jy-state  conditions  were  obtained  by  the  present  authors  [l]. 
By  choosing  the  time  t  =  0  as  the  instant  when  the  radiation  flux  qQ  is  cut 
off,  the  problem  is  reduced  to  finding  the  temperature  and  stress  solutions 
for  a  plate  with  known  initial  temperature  distribution  cooled  in  the  rear 
surface  by  Newtonian  convection  with  heat  transfer  coefficient  h.  The  optical 
properties  of  the  material  are  assumed  to  be  "gray",  i.e.,  independent  of 
wavelength.  The  surface  reflectivity  and  thickness  of  the  plate  are  assumed 
to  be  such  that  the  effect  of  multiple  internal  reflections  can  be  neglected. 

The  material  properties  that  affect  the  temperature  distribution  and  thermal 
stresses,  such  as  the  emissivity  (absorptivity),  absorption  coefficient, 
thermal  conductivity  and  thermal  diffusivity,  coefficient  of  thermal  expansion. 
Young's  modulus  and  Poisson's  ratio  are  assumed  to  be  independent  of  temperature. 


It  is  also  assumed  that  the  temperature  of  the  plate,  over  the  time  period 
required  for  stresses  to  reach  a  value  nearly  equal  to  their  maximum,  remains 
sufficiently  low  so  that  the  reemitted  radiation  is  negligible  in  comparison 
to  the  incident  radiation. 

Solutions  for  the  transient  temperature  T  are  obtained  [3]  by  solving 


3T  =  13T 
3x2  <3t 


where  k  is  the  thermal  diffusivity. 

The  initial  condition  for  the  temperature  distribution  for  the  present 
problem  is  [l], 

00  ^ 

T(x,0)  =  Tq  +  - y  cos  l>n(a+x)]  (2a) 

icAn 

where  the  (An)  are  the  roots  of  the  transcendental  equation 


Antan(2Xna)  =  k 


<q  ep  [p+e  2ya{A  sin(2A  a)-p  cos(2A  a)} 


Gn  N  k 
n 


2  2 
p  +A 

n 


N  =  ttt — [sin(2A  a)cos(2A  a)  +  2A„a] 
n  2A  a  n  n  '  n  J 


The  boundary  conditions  are: 


3T  (  ^ 

3x 


3T.  ^ 

3^(a»t) 


^T(a,t)  -  Tq] 
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Transient  Temperatures 


The  solution  for  the  temperature  was  obtained  by  solving  the  differential 
equation  (1)  subject  to  the  conditions  given  by  Eqs.  (2)  -  (4)  and  is  fonnd 
to  be 


T(x,t)  =  Tq  +  Z=0  Bn(t)  cos  Xn(x+a) 


where 

G  ,  2 

d  /_■>  n  -kA  t 

Bn(t)  =  - 2  6  n 

kA 

n 

with  An  and  Gn  given  by  Eqs.  (2b)  and  (2c), 


Thermal  Stresses 


(5) 


(6) 


Thermal  stresses  are  obtained  by  substitution  for  T  from  Eq.  (5) 
in  [4]. 


aE 

7  ~  i - 

y,z  1-v 


[-T  +  i  / 


Tdx  + 


-a 


V 

2aJ  -a 


Txdx] 
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which  results  in 
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5!  1=0  Bn(t)  “s  [X„(,I+a)  ^ 
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aE 

2a(l-v) 


B„(t) 


sin(2X  a) 
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+  3xaE  f  R  ft.\ 
2a(l-v)  n=0  n(t) 


sin(2A  a) 
n 


aA 
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2 S"  “S  (21na)  '  XT 
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For  convenience,  the  numberical  results  are  reported  in  terms  of  the 
non-dimens ional  tempera ture , 

*  (T-T)k 

T  =  - - — 

eq  a 
o 

the  non-dimensional  stress, 

*  -  s(l~~v)k 

°  aEq^Ea 

A  2 

the  non-dimensional  time  t  =  <t/a  ,  and  the  non-dimensional  distance  x/a. 

NUMBERICAL  RESULTS  AND  DISCUSSION 

The  analytical  results  presented  above  can  be  illustrated  by  a  numerical 
example.  Because  of  the  relatively  large  number  of  independent  variables, 
the  numerical  results  to  be  presented  will  be  limited  to  indicate  a  number 
of  general  results.  The  thermal  conductivity  of  the  plate  was  taken  as  K  = 

0.3  W.cm  1.°C~1,  a  typical  value  for  a  dielectric  material  near  room  tempera¬ 
ture.  Also,  the  half-thickness  of  the  plate,  a,  was  taken  as  1  cm,  except 
where  noted. 

Figure  la  and  lb  show  the  distributions  of  the  temperatures  and  thermal 

stresses  at  a  time  immediately  prior  to  the  interruption  of  radiation  at  t  *  0 

2q  ^ 

for  a  value  of  the  convective  heat  transfer  coefficient  h  =  0.006  W.cm.  C 
These  temperatures  and  stresses,  which  represent  the  initial  conditions  for 
the  present  study,  resulted  from  the  combined  effects  of  radiation  heating 
and  convective  cooling  of  sufficiently  long  duration  that  the  temperatures 
and  stresses  have  reached  steady-state  conditions.  For  many  of  the  details, 
the  reader  is  referred  to  the  earlier  study  of  the  temperature  and  thermal 
stresses  [l]  during  the  heat-up  stage  of  the  plate.  The  steady-state  tempera¬ 
ture  at  the  rear-face  is  an  inverse  function  of  the  heat  transfer  coefficient. 
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The  temperature  profile  and  the  resulting  thermal  stresses  in  the  plate  were 
found  to  be  independent  of  the  value  of  the  heat  transfer  coefficient.  The 
maximum  value  of  thermal  stress  was  found  to  occur  at  a  value  of  optical 
thickness  ya  -  2.  At  ya  =  0  the  stresses  are  zero  because  no  heat  is  absorb¬ 
ed.  As  ya  -*•  ®,  the  temperature  distribution  becomes  linear  which  also 
results  in  a  state  of  zero  thermal  stress.  Due  to  the  downward  curvature 
of  the  temperature  distribution,  the  thermal  stresses  are  tensile  at  the 
hottest  section  of  the  plate  (x  =  -a)  [5], 

Figures  2a  and  2b,  for  a  value  of  optical  thickness  ya  =  3,  show  the 
transient  temperatures  and  stresses  following  the  interruption  of  the  radiation 
at  t  *  0  for  a  number  of  values  of  time,  t,  for  the  initial  conditions  shown 
in  Fig.  la  and  lb.  The  stresses  initially  rise  to  a  maximum  followed  by  a 
decrease  with  time.  The  position  of  the  maximum  compressive  stress  within 
the  plate  shifts  with  increasing  time. 

The  transient  thermal  stresses  at  both  the  front  and  rear  faces  of  the 
plate  are  shown  in  Figs.  3a  and  3b,  respectively,  for  a  number  of  values  of 
the  optical  thickness  for  the  initial  conditions  shown  in  Figs,  la  and  lb. 

It  was  found  that  for  ya  <  0.5,  the  stresses  decreased  with  increasing  time 
without  passing  through  a  maximum. 

Figure  4  presents  evidence  for  the  coupled  nature  of  the  solutions  for  the 
temperature  and  thermal  stress  due  to  the  interruption  of  the  radiation  and  those 
for  the  heating  stage  of  the  plate  prior  to  the  radiation  cut-off.  For  ya  =  3, 
Fig.  4  shows  the  transient  thermal  stresses  at  the  rear  surface  (x  =  a  =  1  cm)  for 
a  range  of  values  of  the  heat  transfer  coefficient.  Even  for  a  range  in  h  of  a 
factor  of  15,  the  maximum  values  of  the  thermal  stresses  differ  by  only  a  few  per¬ 
cent.  At  least  qualitatively,  this  effect  arises  because  the  temperature  rise  of 
the  plate  above  ambient  at  t  -  0  is  inversely  proportional  to  the  heat  trans¬ 
fer  coefficient.  On  subsequent  cooling,  the  magnitude  of  the  thermal  stresses, 
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at  least  for  small  values  of  the  Biot  number  3  ■  ah/k,  are  directly 
proportional  to  the  heat  transfer  coefficient  as  well  as  the  total 
temperature  range  over  which  the  plate  is  being  cooled.  This  leads  to 
the  net  effect  that  the  peak  stresses  following  the  cut-off  radiation, 
to  a  first  approximation,  are  independent  of  the  heat  transfer  coefficient. 

The  stresses  in  the  front  and  rear  surface  for  a  plate  half-thickness 
a  =  10  cm  and  pa  =  3  are  shown  in  Figs.  5a  and  5b,  respectively.  Comparison 
of  the  data  of  Figs.  5a  and  5b  and  Figs.  3a  and  3b  shows  that  for  plates  of 
different  thickness  with  identical  values  of  optical  thickness,  non-dimen- 
sionalization  of  the  numerical  data  does  not  result  in  a  universal  plot 
from  which  the  thermal  stresses  can  be  obtained  for  any  dimension.  This  can 
be  explained  by  noting  that  for  the  same  value  of  optical  thickness  pa,  the 
total  amount  of  heat  absorbed  at  steady-state  (t*0)  is  independent  of  the 
plate  thickness.  For  this  reason  the  temperature  gradient  within  the  plate 
at  the  rear  surface  is  the  same  for  all  thicknesses.  The  temperature  rise 
in  the  front  of  the  plate,  however,  is  not  proportional  to  the  plate  thick¬ 
ness  (unless  pa  -*■  “) .  For  this  reason  the  steady-state  temperature  dis¬ 
tribution  at  t*=  0  does  not  become  normalized  on  a  dimensionless  basis.  This 
can  easily  be  seen  by  re-deriving  Eq.  (2a)  in  another  form;  i.e.,  if  we 
solve  the  steady  heat  conduction  equation  with  the  same  heat  source  used  to 
derive  Eq.  (2a)  Instead  of  taking  the  limit  of  the  transient  solution,  we  find 


[T(x)-TQ]k 


e-2pa_e-px 


♦  +  ±  «  -  *'2ua> 


(9) 


in  non-dimensional  form.  Clearly,  if  a  is  varied  while  pa  is  held  constant, 
the  non-dimensional  temperature  difference  cannot  remain  constant. 
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A  similar  effect  Is  noted  In  Part  II  of  this  study  for  front  face  cooling. 
This  suggests  that  for  plates  of  different  thickness  separate  calculations  of 
the  temperatures  and  stresses  will  be  required. 
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Captions 


Fig.  1.  Distribution  of  steady  state  temperature  (a)  and  thermal  stress  (b) 
in  a  partially  absorbing  flat  plate  heated  in  the  front  by  normally 
incident  radiation  and  cooled  by  convection  at  the  rear  surface  for 
various  values  of  pa  with  k  =  0.3  W.cm.-l°C.~l,  h  =  0.006  W.cm.~2°C.~l, 
k  =  0.1  cm.2s.^l  and  a  =  1  cm. 


Fig.  2.  Distribution  of  temperature  (a)  and  thermal  stress  (b)  for  various 
values  of  time  following  the  interruption  of  the  radiation  at 
time  t  =  0  in  a  partially  absorbing  flat  plate  with  originally  steady 
state  temperature  and  stress  distributions  due  to  normally  incident 
radiation  in  the  front  face  and  cooling  at  the  rear  surface  by  con¬ 
vection  for  k  =  0.3  W.cm.”l°c._l,  h  ■  0.006  W.cm.“2oc.~l,  k  *  0.1 
cm.^s.-l,  pa  =  3  and  a  =  1  cm. 


Fig.  3.  Transient  thermal  stresses  in  the  front  (a)  and  the  rear  face  (b) 
following  the  interruption  of  the  radiation  at  time  t  =  0  in  a 
partially  absorbing  flat  plate  with  originally  steady  state  tem¬ 
perature  and  stress  distribution  due  to  normally  incident  radiation 
in  the  front  face  and  cooling  by  convection  at  the  rear  surface 
for  various  values  of  pa  with  k  *  0.3  W.cm.-loC.-l,  h  *  0.006  W.cm.“2oc.~l, 
k  *  0.1  cm.2s.”l  and  a  *  1  cm. 


Fig.  4.  Transient  thermal  stress  in  the  rear  surface  following  the  interruption 
of  the  radiation  at  time  t  =  0  in  a  partially  absorbing  flat  plate  with 
originally  steady  state  temperature  and  stress  distribution  due  to 
normally  incident  radiation  in  the  front  face  and  cooling  by  convection 
at  the  rear  surface  for  various  values  of  n  with  k  =  0.3  W.cm.~l°C.~l, 
pa  -  3,  kc  *  0.1  cm.2s_1  and  a  =  1  cm. 


Fig.  5.  Transient  thermal  stresses  in  the  front  (a)  and  the  rear  face  (b) 
following  the  interruption  of  the  radiation  at  time  t  =  0  in  a 
partially  absorbing  flat  plate  with  originally  steady  state  temperature 
and  stress  distribution  due  to  a  normally  incident  radiation  in  the 
front  face  and  cooling  by  convection  at  the  rear  surface  for  various 
values  of  pa  with  k  *  0.3  W.cm."l°C.~l,  h  ■  0.006  W.cm.~2oc.-l,  *  = 

0.1  crn.^s.-!  and  a  ■  10  cms. 
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Fig.  1.  Distribution  of  study  stmts  tampsraturs  (m)  mnd  thsrmal 
strsss  (b)  in  m  pmrtimlly  absorbing  flat  pUts  hutsd 
in  ths  front  by  normally  incidsnt  radiation  and  coolsd 
by  convsction  at  ths  rur  surfacs  for  various  valuss  of 
pa  with  k  »  0.3  h  *  0.006  W.ca.”’  °C.”  , 

<  •  0.1  cm^s.*1  and  a  •  1  cm. 
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Fig.  A.  Transient  thermal  stress  in  the  rear  surface  following  the 

Interruption  of  the  radiation  at  time  t  ■  0  in  a  partially  ab¬ 
sorbing  flat  plate  with  originally  steady  state  temperature  and 
stress  distribution  due  to  normally  incident  radiation  in  the 
front  face  and  cooling  by  convection  at  the  rear  surface  for 
various  values  of  h  with  k  -  0.3  W.cm."^°C.“1,  ua  -  3,  <  - 
0.1  cm.^s"!  and  a  •  1  cm. 


Pig.  2.  Distribution  of  temperature  (a)  and  thermal  stress 
(b)  for  various  values  of  time  following  the  inter¬ 
ruption  of  the  radiation  at  time  t  =  0  in  a  partial¬ 
ly  absorbing  flat  plate  with  originally  steady  state 
temperature  and  stress  distributions  due  to  normally 
incident  radiation  in  the  front  face  and  cooling  at 
the  rear  surface  by  convection  for  k  =  0.3  W.cm._-*-°C.-^ , 
h  *  0.006  W. cm. ~^°C . “1 ,  <  =  0.1  cm.^s."!,  pa  =  3  and 
a  ■  1  cm. 
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Transient  thermal  stresses  in  the  front  (a)  and  the  rear  face  (b) 

following  the  interruption  of  the  radiation  at  time  t  =  0  in  a 

partially  absorbing  flat  plate  with  originally  steady  state  tem¬ 
perature  and  stress  distribution  due  to  a  normally  incident  radia 
tion  in  the  front  face  and  cooling  by  convection  at  the  rear  sur¬ 
face  for  various  values  of  pa  with  k  =  0.2  W.cm.'l0C.-1 •  h  =  0.00 

W.cm.“2oc,~l  <  ~  0.1  cm.2s.~l  and  a  *  10  cms. 
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ABSTRACT 


In  Part  II  of  this  study  the  thermal  stresses  are  analyzed  in  a 
partially  absorbing  flat  plate,  due  to  a  sudden  interruption  of  steady- 
state  asymmetric  radiation  on  the  convectively  cooled  front  face. 

As  in  the  case  of  Part  I,  the  maximum  tensile  stresses  occur  in  the  rear 
face.  For  the  optical  thickness  pa  >  0.5,  the  maximum  tensile  stress 
initially  rises  to  a  peak  and  then  decreases  with  time.  On  the  other 
hand,  for  pa  <  0.5,  the  stresses  decrease  monotonically  with  time. 

A  comparison  of  the  results  for  the  front  (Part  II)  and  rear  (Part 
I)  cooling  indicates  that  values  of  the  maximum  stress  are  numerically 
identical  for  small  values  of  heat  transfer  coefficient,  h,  and  plate 
thickness,  a.  Increasing  difference  in  the  magnitude  of  the  stresses 
for  the  two  cases  is  found  with  increasing  h  and  a. 


INTRODUCTION 


In  Part  I  of  this  study  [1],  solutions  were  presented  for  the  tempera¬ 
tures  and  thermal  stresses  due  to  the  sudden  interruption  of  long-term 
steady-state  asymmetric  thermal  radiation  on  the  front  surface  of  a  partially 
absorbing  flat  plate  cooled  by  Newtonian  convection  at  the  rear  surface. 

In  part  II  the  same  problem  is  considered,  with  the  exception  that  the 
radiation  heating  and  convection  cooling  both  occur  at  the  same  (front) 
surface  of  the  plate.  The  assumptions  and  notation  for  the  present  problem 
are  identical  to  those  of  Part  I. 
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Solutions  for  the  transient  temperature  distribution  are  ohtalned 
by  solving: 


32T  ^  1  9T 

9x2  "  *  3t 

subject  to  the  boundary  conditions: 


3T  ,  . 

to  (a’fc)  "  0 

to  (‘a’t)=~k  fT(-a,t)-To] 


(2) 

(3) 

(4) 


Transient  Temperatures 

Solution  of  Eq .  2  subject  to  the  initial  and  boundary  conditions 
given  by  Eqs.  (1,3,4)  yields: 


T(x,t)  =  Tq  +  Bn(t)  cos  Xn  (a-x) 


where 
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with  X  and  G  defined  earlier, 
n  n 
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Thermal  Stresses 


Following  the  procedure  outlined  in  Part  I  the  thermal  stresses 

are: 
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The  numerical  results  will  be  presented  in  terns  of  the  non-dimen¬ 
sional  temperature,  stress  and  tine  defined  in  Part  1. 

NUMERICAL  RESULTS  AND  DISCUSSION 

In  order  to  facilitate  comparison,  the  data  presented  will  be  for 
the  identical  values  of  the  pertinent  parameters  chosen  for  Part  I,  namely 
k  *  0.3  W.cm.  *°C  *  and  a  -  1  cm,  except  where  noted.  Furthermore,  the 
sequence  of  the  figures  will  be  identical.  For  a  range  of  values  of 
optical  thickness,  pa,  Figs,  la  and  lb  show  the  distribution 
of  the  temperature  and  thermal  stresses,  respectively,  in  the 
plate  at  t  *  0,  immediately  prior  to  the  interruption  of  the  radiation. 
Comparison  of  the  temperature  distribution  (Fig.  la)  with  the  corresponding 
data  in  Part  I  shows  that  radiation  heating  and  convective  cooling  at  the 
same  surface  produces  lower  temperatures  at  the  exposed  surface  than 
for  radiation  heating  and  convective  cooling  at  opposite  surfaces.  As 
pointed  out  in  earlier  studies  12,3],  this  latter  conclusion  could  be 
relevant  to  the  design  of  components  of  high-intensity  solar  collectors 
subject  to  material  degradation.  Excessive  surface  temperatures  can  he 
avoided  by  providing  the  convective  cooling  at  the  same  surface  being 
irradiated. 

The  values  of  the  thermal  stresses  reach  their  maximum  for  a  value  of 
optical  thickness  pa  -  2.  Comparison  of  the  values  of  the  thermal  stresses 
at  t  ■  0,  just  prior  to  the  interruption  of  the  radiation  as  displayed  in 
Fig.  lb  of  parts  I  and  II,  shows  them  to  be  almost  independent  of  the  surface 
at  which  the  convective  cooling  takes  place.  However,  comparison  of  the 
analytical  expressions  for  the  thermal  stresses  for  the  two  modes  of 
cooling  show  them  not  to  be  the  same.  Further  detailed  analysis  reveals 


that  the  numerical  values  for  the  thermal  stresses  for  given  values  of 
the  heat  transfer  coefficient,  optical  thickness,  and  thermal  conductivity 
become  equal  only  when  the  plate  thickness  2a  0. 

Figures  2a  and  2b,  for  pa  *  3,  show  the  distribution  of  the  transient 
temperatures  and  stresses,  respectively  for  a  range  of  values  of  time,  x, 
following  the  interruption  of  the  radiation  using  the  same  values  of  k,  k  and 
h  as  in  Figs,  la  and  lb.  The  data  for  the  stress  are  almost  identical 
to  those  given  in  Part  I.  However,  slight  differences  necessitated  choosing 
different  values  for  t  in  order  to  avoid  too  much  overlap  of  the 
curves,  for  purposes  of  clarity. 

Figures  3a  and  3b  show  the  values  for  the  thermal  stresses  at  the 
front  and  rear  faces  of  the  plate,  respectively,  as  a  function  of  time 
for  a  range  of  values  of  pa.  As  found  in  Part  I,  the  stresses  for 
pa  ~  0.5  initially  increase  to  a  maximum  followed  by  a  monotonic  decrease 
with  increasing  time.  The  stresses  in  the  rear  face  exceed  those  in  the 
front  face. 

In  Fig.  4,  we  compare  the  transient  thermal  stresses  at  the  rear  face 
for  a  range  of  values  of  the  heat  transfer  coefficient  for  pa  =  3.  At 
thermal  equilibrium,  just  prior  to  the  cut-off  of  the  radiation,  the 
magnitude  of  stress  is  independent  of  the  heat  transfer  coefficient  [2], 

This  effect  arises  because  the  heat  transfer  coefficient  controls  the 
absolute  temperature  at  the  front  face,  but  not  the  internal  temperature 
distribution,  which  is  responsible  for  the  development  of  the  thermal 
stresses.  With  increasing  time,  the  stresses  go  through  a  maximum,  with 
a  peak  value  which  is  an  inverse  function  of  but  strongly  dependent  on 
the  magnitude  of  the  heat  transfer  coefficient. 


Comparison  of  the  present  data  in  Fig.  4  with  the  corresponding  data 

for  rear  face  cooling  presented  in  Part  I  shows  that  for  h  =  0.006  cal. 

-2  -1  -1 

cm.  °C.  s  ,  the  numerical  values  for  the  stresses  at  any  instant  of 

time  are  almost  identical.  In  contrast,  at  the  highest  value  of  the 

-2  -1  -1 

heat  transfer  coefficient  (h  *  0.09  cal. cm.  °C.  s  ),  considerable 
difference  exists  in  the  magnitude  of  the  stresses  for  front  and  rear 
face  cooling.  This  difference  was  sufficiently  large  that  a  change  in 
scale  was  required  for  presentation  of  the  respective  data  shown  in  Fig.  4 
in  Parts  I  and  II. 

The  difference  in  the  effect  of  the  heat  transfer  coefficient  on 
the  magnitude  of  thermal  stress  for  front  and  rear  face  cooling  can  be 
explained  qualitatively.  It  should  be  noted  that  the  magnitude  of  the 
transient  thermal  stresses  for  both  front  and  rear  face  cooling  will 
be  a  function  of  both  the  magnitude  of  the  initial  temperature  rise 
of  the  plate  as  well  as  the  temperature  non-uniformity  within  the 
plate.  As  noted  earlier  the  initial  temperature  rise  of  the  plate  is 
inversely  proportional  to  the  value  of  the  heat  transfer  coefficient. 

The  non-uniformity  of  the  internal  temperature  distribution  is  independent 
of  the  value  of  the  heat  transfer  coefficient  and  is  a  function  of  the 
distribution  of  the  internal  absorption  of  the  incident  radiation.  Two 
limiting  cases  can  be  considered,  namely  h  •*  0  and  h  -*■  <*>,  For  h  •+•  0, 
the  temperature  non-uniformity  within  the  plate  will  become  very  small 
relative  to  the  initial  total  temperature  rise.  This  latter  quantity  for 
front  and  rear  face  cooling  will  tend  to  become  equal.  With  the  magnitude 
of  the  transient  thermal  stresses  being  proportional  to  this  initial 
temperature  rise,  the  stresses  will  be  equal  for  front  and  rear  face 
cooling. 


As  h  -*■  »,  a  completely  different  situation  arises.  For  front  face 
cooling,  no  temperature  rise  and  temperature  non-uniformity  occurs 
which  yields  zero  initial  and  transient  thermal  stress  on  cooling.  Rear 
face  cooling,  however,  will  result  in  a  temperature  rise  at  the  front 
face  and  accompanying  non-uniform  temperature  distribution  within  the  plate. 
In  transient  cooling  after  the  cut-off  of  the  radiation,  transient  thermal 
stresses  will  inevitably  result.  Obviously,  the  results  for  the 
present  study  (0  <  h  <  ®)  will  fall  between  the  above  mentioned  two 
limits. 

Figures  5a  and  5b  show  the  transient  thermal  stresses  at  the  front 
and  rear  surface  of  the  plate,  respectively,  for  a  =  10  cm  for  the  same 
numerical  values  of  the  appropriate  parameters  for  the  data  in  Figs. 

3a  and  3b.  Comparison  with  these  latter  data  indicates  that  (as 
also  noted  in  Part  I)  plotting  the  data  for  the  present  thermal 
problem  on  a  dimensionless  basis  does  not  result  in  a  universal  plot 
for  the  same  reasons  as  outlined  in  Part  I. 

Comparison  of  Figs.  5a  and  5b  in  Part  1  and  II  shows  that  for  a  =  10  cm, 
the  peak  transient  thermal  stresses  for  front  face  cooling  lie  somewhat 
below  the  corresponding  values  for  rear  face  cooling.  As  noted  earlier, 
for  a  *  1  cm,  these  stresses  are  almost  identical.  It  appears  that  the 
magnitudes  of  the  peak  tensile  stresses  for  front  and  rear  face  cooling 
diverge  with  increasing  plate  thickness. 

In  summary,  solutions  and  numerical  values  have  been  presented  for 
the  transient  thermal  stresses  in  a  semi-absorbing  flat  plate  initially 
at  thermal  equilibrium  due  to  long-term  radiation  heating  and  convective 
cooling  at  the  front  surface,  and  subjected  to  an  instantaneous  cut-off 


of  the  radiation.  The  transient  stresses  were  found  to  be  quantitatively 
almost  identical  to  those  for  convective  cooling  at  the  rear  surface, 
for  small  plate  thicknesses  and  low  heat  transfer  coefficient.  However, 
the  magnitude  of  the  transient  stresses  for  front  and  rear  cooling 
diverge  with  increasing  plate  thickness  and  heat  transfer  coefficient. 
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CAPTIONS 


Fig.  I. 


Fig.  2. 


Fig.  3. 


Fig.  4. 


Distribution  of  steady  state  temperature  (a)  and  thermal 
stress  (b)  in  a  partially  absorbing  flat  plate  heated  by 
normally  incident  radiation  and  cooled  by  convection  in 
the  front  face  for  various  values  of  pa  with  k  =  0.3 
W.cm.-l°C"l ,  h  =  0.006  W.cm.~2°C."^,  k  =  0.1  cm.s-^  and 
a  =  1  cm. 

Distribution  of  temperature  (a)  and  thermal  stress  (b) 

for  various  values  of  time  folloiwng  the  interruption  of 

the  radiation  at  time  t  =  0  in  a  partially  absorbing  flat 

plate  with  originally  steady  state  temperature  and  stress 

distributions  due  to  normally  incident  radiation  heating 

and  convection  cooling  at  the  front  face  for  k  ■  0.3  w.cm."  °C  , 

h  =  O.OOfi  W. cm.~^0C-^,  k  =  0.1  cm.^s.-^,  pa  =  3  and  a  =  1  cm. 

Transient  thermal  stresses  in  the  front  (a)  and  the  rear 
face  (b)  following  the  interruption  of  the  radiation  at 
time  t  =  0  in  a  partially  absorbing  flat  plate  with 
originally  steady  state  temperature  and  stress  distribution 
due  to  normally  incident  radiation  heating  and  convection 
cooling  at  the  front  face  for  k  =  0.3  W.cm.-l°C.-l ,  h  = 

0.006  W.cm.“2°C.-l,  <  *  0.1  cra.^s.-!  and  a  =  1  cm. 


Transient  thermal  stresses  in  the  rear  surface  following 
the  interruption  of  the  radiation  at  time  t  =  0  in  a  partially 
absorbing  flat  plate  with  originally  steady  state  temperature 
and  stress  distribution  due  to  normally  incident  radiation 


heating  and  convection  cooling  in  the  front  face  for  various 
values  of  h  with  k  =  0.3  W.cm.~l°C.~l,  pa  =  3,  <  =  0.1 
cm.^s.-l  and  a  =  1  cm. 


Fig.  5.  Transient  thermal  stresses  in  the  front  (a)  and  the  rear  face 
(b)  following  the  interruption  of  the  radiation  at  time  t  =  0 
in  a  partially  absorbing  flat  plate  with  originally  steady- 
state  temperature  and  stress  distribution  due  to  a  normally 
incident  radiation  heating  and  convection  cooling  in  the 
front  face  for  various  values  of  pa  with  k  =  0.3  W.cm."^°C.~^ , 
h  =  0.006  W.cm.-20c-l,  <  =  0.1  cm.^s.-!  and  a  =  10  cms. 


NON -DIMENSIONAL  DISTANCE,  n/o 
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Fig.  2.  Distribution  of  temperature  (a)  and  thermal  stress  (b) 
for  various  values  of  time  following  the  interruption 
of  the  radiation  at  time  t  =  0  in  a  partially  absorbing 
flat  plate  with  originally  steady  state  temperature  and 
stress  distribution  due  to  normally  incident  radiation 
heating  and  convection  cooling  at  the  front  face  for 
k  =  0.3  W.  cm.  -loC.~l ,  h  =  0.006  W.cm.-^C'1 ,  k-  =  0,1  cm.2 
s."l,  pa  =  3  and  a  =  1  cm. 


Fig.  3.  Transient  thermal  stresses  in  the  front  (a)  and  the  rear  face 
(b)  following  the  interruption  of  the  radiation  at  time  t  =  0 
in  a  partially  absorbing  flat  plate  with  originally  steady 
state  temperature  and  stress  distribution  due  to  normally 
incident  radiation  heating  and  convection  cooling  at  the 
front  face  for  k  =  0.3  W.cm.-1°C-1,  h  =  0.006  W.cm.-2°C.-1 , 
k  =  0.1  cm.2s.“l  and  a  =  1  cm. 
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Transient  thermal  stress  in  the  rear  surface  following  the 
interruption  of  the  radiation  at  time  t  =  0  in  a  partially 
absorbing  flat  plate  with  orginally  steady  state  tempera¬ 
ture  and  stress  distribution  due  to  normally  incident 
radiation  heating  and  convection  cooling  in  the  front  face 
for  various  values  of  h  with  k  =  0.3  W.cm.“l°C""l,  ya  =  3, 
k  =  0.1  cm.2s.~l  and  a  =  1  cm. 
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Fig.  5.  Transient  thermal  stresses  in  the  front  (a)  and  the  rear  face 
(hi  following  the  interruption  of  the  radiation  at  time  t  =  0 
in  a  partially  absorbing  flat  plate  with  originally  steady- 
state  temperature  and  stress  distribution  due  to  a  normally 
incident  radiation  heating  and  convection  cooling  in  the  front 
face  for  various  values  of  pa  with  k  =  0.3  W.cm.-1°C-1 ,  h  = 
0.006  W.cm.  ^  c~l ,  k  *  0.1  cm.^g.-l  and  a  =  10  cms. 
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In  a  recent  study  (1),  the  present  authors  presented  an  analysis  of  the 
frequency  dependence  of  the  thermal  stresses  in  a  partially  absorbing  flat 
plate  due  to  the  internal  absorption  of  cyclic  thermal  radiation  normally 
incident  on  one  face,  and  cooling  by  convection  at  the  opposite  face.  At 
higher  frequencies,  after  an  initial  transient,  the  plate  reaches  a  steady- 
state  temperature  profile  and  associated  thermal  stress  distribution  with  a 
superposed  cyclic  variation  of  small  amplitude  relative  to  the  total  temperature 
rise  and  thermal  stresses. 

At  very  low  frequencies,  a  different  effect  occurs.  The  mean  temperature 
of  the  plate  and  associated  thermal  stresses  cycle  along  with  the  intensity  of 
the  incident  heat  flux.  As  the  frequency  w  -*o,  the  plate  is  subjected  to  a 
quasi  steady-state  heat  flux.  This  effect  could  not  be  included  in  the  linear 
plots  of  the  frequency  dependence  of  the  maximum  thermal  stresses  encountered 
during  the  first  cycle  and  after  many  cycles  at  which  the  plate  has  reached  a 
constant  time-averaged  temperature  given  in  Figs.  5a  and  5b,  respectively,  of 
the  original  study  (1).  The  purpose  of  this  communication  is  to  report  these 
values.  Material  property  values  identical  to  those  used  for  the  data  given 
in  Fig.  5  of  the  original  paper  were  used  for  the  calculation,  namely:  heat 
transfer  coefficient  h  =  0.006  W.cm^,°C"\  thermal  diffusivity  <  =  0.1  cm^.s"\ 


plate  hal f- thickness  a  =  1  cm  and  absorption  coefficient  u  =  3  cm"^. 

Figures  la  and  lb  show  the  values  of  the  thermal  stresses  in  the  front 
face,  center  and  back  face(x  =  -a,  0,  and  a,  respectively)  as  a  function  of 
the  angular  frequency  encountered  during  the  first  cycle  and  after  the  plate 
has  reached  a  quasi-static  equilibrium  when  no  cycle-to-cycle  variation  in 
temperature  or  stress  is  found.  The  results  are  presented  in  terms  of  the 
maximum  and  minimum  stress  encountered  within  a  given  cycle  for  the  three 
positions  in  the  plate. 

Figure  la  shows  that  as  w  -*  o,  the  mimimum  value  of  tensile  stress  of 
0.334  encountered  during  the  first  cycle  is  twice  the  value  of  0.167  after  long¬ 
term  steady-state  radiation  of  intensity  qQ,  found  in  an  earlier  study  (2).  It 
should  be  noted  that  for  the  specific  cyclic  variation  chosen  for  this  study  with 
q  =  qQ(l  +  sin  wt),  the  maximum  intensity  of  radiation  within  any  one  cycle 
reaches  a  value  of  2qQ.  This  indicates  that  indeed  as  u  +  o  the  plate  can  be 
considered  to  be  subjected  to  a  quasi -steady-state  heat  flux  after  an  initial 
brief  transient  at  t  =  0  when  the  heat  flux  changes  discontinuously  from  zero  to 
qQ.  The  stresses  due  to  this  initial  transient  are  identical  to  those  found 
for  sudden  imposition  of  steady-state  radiation  with  a  constant  flux  qQ  (2). 

Figure  lb  shows  that  after  long-term  radiation,  as  u>  o  the  stresses  in 
all  three  positions  can  become  zero.  This  is  expected  since  within  the  cycle 
the  value  of  minimum  heat  flux  is  identically  equal  to  zero.  The  other  (non¬ 
zero)  values  of  stress,  as  u  +  o,  corresponds  to  the  maximum  value  of  heat  flux 
within  the  cycle,  equal  to  2qQ. 

Of  interest  to  note  is  that,  as  indicated  by  both  figures,  the  magnitude 
of  the  maximum  stress,  as  well  as  the  position  within  the  plate  at  which  it  occurs, 
varies  with  frequency.  When  the  frequency  «o  -*■  ®  the  stresses  correspond  to 
those  found  for  a  constant  heat  flux  qQ. 
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1,  Frequency  dependence  of  range  of  thermal  stresses  in 
front,  center  and  rear  in  partially  absorbing  flat 
plate  heated  in  front  by  cyclic  thermal  radiation  and 
cooled  by  convection  at  the  rear  surface  for:  h  = 
0.006  W.cm.-2°C'l,  k  =  0.3  W.cm.-'OC"1,  <  =  0.1  cm.^s 
y  =  3  cm"1  and  a  =  1  cm.  (a)  during  first  cycle  and 
(b)  for  steady  state. 


CHAPTER  IX 


Analysis  of  Thermal  Stress  Failure 
of  Segmented  Thick-Walled  Refractory  Structures 

by 

W.  S.  Chang,  C.  E.  Knight,  D.  P.  H.  Hasselman,  and  R.  G.  Mitchiner 


Departments  of  Mechanical  and  Materials  Engineering 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061  USA 


ABSTRACT 


A  finite-element-analysis  was  conducted  of  the  thermal  stresses  in 
rectangularly  shaped  components  of  segmented  thick-walled  refractory 
structures  with  length  in  the  direction  of  heat  flow  much  larger  than 
the  height  or  width.  Heating  was  assumed  to  occur  uniformly  at  a  linear 
rate  over  the  front  face  defined  by  the  smaller  dimensions.  Most  cal¬ 
culations  were  carried  out  for  two-dimensional  conditions  of  plane 
stress,  which  assumes  that  the  magnitude  of  stress  is  independent  of 
height.  A  limited  number  of  three-dimensional  calculations  permitted 
analyzing  the  effect  of  height. 

For  moderate  values  of  heating  rates,  corresponding  to  those  pos¬ 
sibly  encountered  in  practice,  a  biaxial  tensile  stress  was  found  to 
exist  some  distance  behind  the  hot  face  with  the  maximum  value  of  stress 
along  the  centerline  parallel  to  the  component  length  in  accordance  with 
the  fracture  mode  observed  in  practice.  At  high  heating  rates,  however, 
the  maximum  tensile  stress  component  was  found  to  occur  parallel  tc  the 
face  being  heated. 

For  conditions  of  plane  stress  under  moderate  heating  rates,  the 
magnitude  of  stress  was  found  to  be  proportional  to  heating  rate  and  an 
inverse  function  of  the  thermal  diffusivity.  Maximum  values  of  stress 
were  encountered  for  intermediate  values  of  width.  This  latter  result 
suggests  that  the  incidence  of  thermal  stress  failure  may  be  reduced  by 
either  reductions  or  increases  in  the  values  of  width  commonly  used  in 
practice.  The  three-dimensional  calculations  indicated  that  the  magni¬ 
tude  of  stress  rises  rapidly  with  height.  This  result  suggests  that 


ii 


improvements  in  thermal  stress  resistance  can  be  obtained  by  reductions 
in  the  values  of  height  from  the  current  values  in  practice. 

By  means  of  a  fracture-mechanical  analysis  it  is  demonstrated  that 
cracks  oriented  perpendicularly  to  the  direction  of  heat  flow  formed  as 
the  result  of  the  above  stresses  could  be  unstable  and  are  expected  to 
propagate  across  the  total  crossection  to  lead  to  total  separation,  in 
accordance  with  observations. 

The  stress  values  obtained  in  the  present  study  were  compared  with 
those  predicted  by  the  Kienow  theory.  Reasonable  agreement  was  found 
for  the  peak  stresses  at  the  lower  values  of  time.  However,  with  in¬ 
creasing  duration  and  the  higher  values  of  height  the  Kienow  theory 
underestimates  the  magnitude  of  stress  considerably. 
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INTRODUCTION 


Engineering  structures  or  components  which  operate  at  elevated  tem¬ 
peratures  require  materials  of  construction  with  high  melting  points, 
combined  with  excellent  resistance  to  corrosion,  deformation  and  frac¬ 
ture  at  high  temperature.  Many  of  these  requirements  are  met  by  engi¬ 
neering  ceramics  or  refractories.  Unfortunately,  these  materials  ex¬ 
hibit  a  high  degree  of  brittleness. 

In  the  majority  of  service  requirements,  involving  elevated  temper¬ 
ature,  the  structures  or  components  inevitably  are  subjected  to  trans¬ 
ient  or  steady-state  heat  flow.  With  few  exceptions  the  corresponding 
spatially  non-uniform  temperature  distribution  within  the  structure  or 
component  can  result  in  the  formation  of  thermal  stresses  of  high  mag¬ 
nitude^.  Because  of  the  high  brittleness  of  most  ceramics  or  refrac¬ 
tories,  failure  under  the  influence  of  these  thermal  stresses  can 
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involve  extensive  crack  propagation  .  Under  these  conditions  failure 
can  be  near  if  not  totally  catastropic  with  an  associated  major  adverse 
effect  on  the  subsequent  satisfactory  performance  of  the  structure  or 
component  in  its  intended  application. 

A  number  of  different  modes  of  failure  and  associated  crack  propaga¬ 
tion  in  brittle  ceramics  and  refractories  can  be  identified.  On  heat¬ 
ing,  corners  of  angularly  shaped  components  can  be  removed  by  spalling 
by  what  appears  to  be  failure  in  shear.  On  heating  or  cooling,  failure 
most  commonly  occurs  in  tension  with  the  direction  of  crack  propagation 

parallel  to  the  direction  of  heat  flow.  This  mode  of  failure  is  entirely 
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in  accordance  with  the  fracture  mode  inferred  from  theoretical  solutions 
for  the  stresses  in  thin  plates,  tubes  and  other  similar  geometries 
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which  indicate  that  the  magnitude  of  the  stresses  in  the  direction  of 
heat  flow  due  to  the  presence  of  the  free  surfaces  are  small  or  even 
negligible  compared  to  the  stresses  perpendicular  to  the  heat  flow.  An 
exception  to  this  general  rule  is  found  in  the  failure  of  thick-walled 
structures  which  consist  of  an  assembly  of  sub-components  with  dimen¬ 
sions  in  the  direction  of  heat  flow  larger  than  those  in  the  perpendicu¬ 
lar  directions.  Examples  of  such  structures  are  linings  of  refractory 
process  vessels  used  in  the  chemical  and  metal  processing  industries, 
coal-gasification  systems,  ship-board  incinerators,  air-preheaters  for 
wind-tunnels,  first  walls  in  fusion  reactors,  and  similar  structures. 

The  linings  of  the  basic-oxygen-furnace  (BOF)  represent  a  particular 
case.  The  dimensions  of  the  individual  components  of  such  linings 
typically  have  cross  sections  measuring  some  8  by  15  cm  with  a  length  of 
as  much  as  75  cm,  with  the  direction  of  heat  flow  along  the  length. 

In  such  individual  components  due  to  the  lack  of  structural  inter¬ 
action  from  component  to  component,  thermal  stresses  perpendicular  to 
the  direction  of  heat  flow  cannot  fully  develop.  In  this  case,  failure 

is  unlikely  to  occur  with  the  plane  of  crack  propagation  parallel  to  the 
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direction  of  heat  flow.  Instead,  as  identified  by  Kienow  and  Ainsworth 
failure  occurs  with  the  plane  of  crack  propagation  oriented  perpendicu¬ 
larly  to  the  direction  of  heat  flow.  In  particular,  this  mode  of  failure 
occurs  during  too  rapid  heat-up  of  the  refractory  lining.  The  origin  of 
failure  appears  to  be  located  on  the  centerline  of  the  refractory  com¬ 
ponent  and  some  5-10  cm  behind  the.  face  at  which  the  heating  occurs. 
Simultaneous  failure  of  a  number  of  adjacent  components  can  result  in 
the  reduction  in  thickness  of  the  refractory  linings  over  a  considerable 


area.  In  turn  this  can  have  a  major  adverse  effect  on  the  total  life¬ 
time  of  the  refractory  structure.  It  should  be  noted  that  a  similar 
mode  of  failure  can  occur  on  cooling^,  with  the  exception  that  the 
origin  of  failure  appears  to  be  located  along  the  side  of  the  component 
with  the  direction  of  crack  propagation  inward  towards  the  center  of  the 
component.  This  study  will  concentrate  on  the  mode  of  failure  during 
heating. 

For  purposes  of  improved  engineering  design  and  assuring  long-term 

structural  reliability  it  is  imperative  that  the  fundamental  variables 

which  control  the  thermal  stress  failure  on  heating  as  described  above 

are  well  understood.  In  an  earlier  study  carried  out  for  this  purpose, 
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Kienow  developed  a  theory  which  relates  the  tensile  stress  in  the 
refractory  component  at  any  position  along  the  center  line  to  the  curva¬ 
ture  (i.e.,  the  second  derivative)  of  the  temperature  distribution  at 
that  position.  Implicit  in  this  relationship  is  that  the  thermal 
stresses  at  any  position  within  the  component  are  independent  of  the 
temperature  distribution  at  other  positions.  This,  however,  is  not 
likely  to  be  correct  in  view  of  the  coupled  nature  of  thermal  stresses 
and  their  general  dependence  on  the  temperature  distribution  and  boundary 

conditions  throughout  a  component  or  structure  as  evidenced  by  an  exten- 
3 

sive  literature  .  In  the  light  of  results  obtained  in  this  study,  it 
should  be  noted  also  that  the  Kienow  theory  predicts  that  the  magnitude 
of  the  stresses  are  a  function  of  the  width  but  not  the  height  of  the 
refractory  component. 

Because  of  the  existence  of  corners  and  free  surfaces  parallel  to 
the  direction  of  heat  flow  analytical  methods  appropriate  for  geometries 
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such  as  infinite  plates,  hollow  and  solid  cylinders  etc.,  are  inappropri¬ 
ate  for  a  theoretical  study  of  the  mode  of  thermal  fracture  of  compo¬ 
nents  of  finite  size,  as  described  above.  In  this  respect,  then,  the 
objective  of  this  study  was  to  investigate  the  above  problem  by  numer¬ 
ical  means,  specifically  the  finite  element  method. 

II.  Numerical  Procedures  and  Boundary  Conditions 

The  mechanical  model  selected  for  this  study,  shown  schematically  in 
Fig.  1,  consisted  of  a  rectangularly  shaped  solid  located  at  0  _<  x  L, 
-b/2  <_  y  £  b/2,  and  -d/2  £  z  _<  d/2  where  L,  b,  and  d  are  the  length, 
width,  and  height,  respectively.  These  dimensions  were  such  that  L  >> 
b,d.  Heating  was  chosen  to  take  place  uniformly  on  the  surface  defined 
by  the  width,  b,  and  height,  d,  located  at  x  =  0.  Because  of  its  signifi¬ 
cance  to  industrial  practice,  this  surface  was  assumed  to  be  heated  at  a 
constant  linear  rate  from  an  initial  uniform  temperature  Tq  at  t  =  0: 

T (x  =  0)  =  at  +  Tq(x  =  0)  (1) 

where  a  is  the  rate  of  change  of  temperature  and  t  is  the  time.  For 
simplicity  Tq  was  taken  equal  to  zero.  The  assumption  of  a  linear  rate 
of  heating  was  also  chosen  for  the  study  of  Ainsworth^.  Heat  flow  was 
assumed  to  be  unidirectional  along  the  x-axis  only.  As  the  segment  is 
assembled  with  other  segments  surrounding  it  there  is  no  heat  transfer 
from  the  sides.  The  rear  surface  at  x  =  L  was  assumed  to  be  in¬ 
sulated,  that  is: 


S  <*  '  L>  ’  0 


(2) 


The  assumption  of  an  insulated  back  surface  is  not  critical  because 
as  found  in  practice,  fracture  occurs  during  the  initial  stages  of 
heating  well  before  the  rear  surface  (x  =  L)  sees  any  perceptible  in¬ 
crease  in  temperature. 

A  major  simplification  in  the  total  computational  effort  was  a- 
chieved  by  assuming  that  the  thermal  conductivity  and  thermal  diffu- 
sivity  were  independent  of  temperature  and  position.  By  this  assump¬ 
tion,  the  transient  temperature  distribution  can  be  obtained  from  an 
analytical  solution  available  from  the  literature,  such  that  only  the 
thermal  stresses  need  to  be  evaluated  numerically. 

As  presented  by  Carslaw  and  Jaeger^  the  transient  temperature  dis¬ 
tribution  for  the  geometry  and  boundary  conditions  described  above  is: 


T(x, t)  -  at  +  +  iSal  £ 
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where 


B  =  x(2n+l) 


(4) 


and  <  is  the  thermal  diffusivity. 

The  thermal  stress  field  in  the  rectangular  solid  is  three-dimen¬ 
sional  even  though  the  temperature  distribution  is  one-dimensional.  The 
stresses  may  be  analyzed  using  three-dimensional  finite  element  models, 
but  for  parametric  studies  such  as  the  current  one,  such  analyses  become 
very  time-consuming  and  expensive.  For  this  reason,  the  stresses  were 

calculated  for  a  two-dimensional  (plane  stress)  distribution  by  means  of 
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a  finite  element  program  used  for  a  number  of  earlier  studies  .  A 
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limited  number  of  three-dimensional  cases  were  also  run  to  investigate 

12 

the  effect  of  component  height,  d,  using  the  SUPERB  finite  element 
program. 

The  two-dimensional  model  consisted  of  a  rectangular  mesh  with  five 
elements  across  the  half  width  by  forty-nine  elements  along  the  length 
and  heat  flow  direction.  The  half  width  model  used  prescribed  boundary 
conditions  on  the  center  line  edge  to  maintain  the  line  of  symmetry. 

For  the  two-dimensional  model,  a  condition  of  plane  stress  were  assumed 
throughout.  This  implies  that  the  results  obtained  are  appropriate  for 
a  plate  with  a  height  much  less  than  its  width. 

Table  1  lists  the  property  values  (or  range  thereof)  used  for  the 
calculations.  The  values  listed  are  typical  of  those  encountered  in 
practice  for  magnesia-based  refractories.  The  influence  of  the  various 
parameters  which  affect  the  magnitude  of  thermal  stresses,  such  as  the 
thermal  diffusivity  or  the  structure  dimensions,  is  not  readily  ex¬ 
pressed  in  dimensionless  form  as  is  common  practice  in  the  analysis  of 
thermal  stresses.  For  this  reason,  the  results  obtained  are  reported  in 
terms  of  absolute  values  of  the  stresses  and  other  pertinent  variables. 
Because  the  magnitude  of  the  stresses  is  directly  proportional  to  the 
values  of  the  coefficient  of  thermal  expansion  and  Young's  modulus  of 
elasticity,  the  stresses  can  be  converted  easily  for  other  values  for 
those  two  properties.  In  order  to  keep  the  total  number  of  stress  cal¬ 
culations  within  reasonable  bounds  the  results  to  be  reported  will 
emphasize  the  general  trend  of  the  role  of  the  various  parameters  on  the 
magnitude  of  thermal  stresses.  Since  brittle  materials  are  most  sus¬ 
ceptible  to  thermal  stress  failure  in  tension,  the  magnitude  of  the 
tensile  thermal  stresses  will  receive  major  emphasis. 
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III.  Results  and  Discussion 


Figs.  2a,  2b,  and  2c  show  typical  distributions  for  the  transient 
temperatures  and  components  of  thermal  stress  in  the  longitudinal  (x) 

and  transverse  (y)  directions,  respectively,  for  numerical  values  of  the 
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thermal  diffusivity,  k  =  12.9  x  10  cm  /s,  heating  rate,  a  *  300°C/hr, 
time,  t  =  3000  s,  width,  b  =  15  cm  and  length,  L  =  60  cm.  For  the  value 
of  time  chosen  for  this  example,  the  back  surface  exhibits  negligible 
increase  in  temperature.  The  surface  being  heated  is  in  a  state  of 
transverse  compression  as  expected.  A  biaxial  tensile  stress  state  is 
found  in  the  interior  some  distance  behind  the  face  being  heated  with 
maximum  values  of  stress  on  the  centerline.  For  this  specific  numerical 
example  the  longitudinal  stress  parallel  to  the  direction  of  heat  flow 
exceeds  the  corresponding  value  perpendicular  to  the  direction  of  heat 
flow.  Such  a  stress  distribution  is  expected  to  result  in  the  mode  of 
failure  observed  in  practice  in  which  the  plane  of  crack  propagation  is 
parallel  to  and  some  distance  behind  the  face  being  heated. 

However,  for  shorter  times  and/or  higher  heating  rates  the  tensile 
stress  state  can  be  such  that  the  transverse  tensile  stresses  exceed  the 
longitudinal  stresses  as  shown  in  Figs.  3a,  3b,  and  3c,  respectively, 
for  a  heating  rate,  a  =  300°C/hr,  time,  t  =  500  s,  and  values  for  the 
other  pertinent  variables  identical  to  those  used  for  the  calculations 
of  the  data  presented  in  Fig.  2.  For  this  reason,  a  tensile  failure 
mode  with  the  plane  of  crack  propagation  parallel  to  the  direction  of 
heat  flow  can  be  anticipated  for  shorter  time  intervals  and  stresses  of 
sufficient  magnitude  as  may  occur  for  the  condition  of  very  rapid  heating 
rates  well  in  excess  of  300°C/hr.  Calculations  of  the  tensile  thermal 
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stresses  in  the  interior  for  other  values  of  width  and  parameters  for 
the  data  shown  in  Figs.  2  and  3  indicated  that  they  always  exhibited 
their  maximum  value  along  the  centerline.  Furthermore,  observed  in 
practice,  failure  always  occurred  parallel  to  the  face  being  heated. 

For  these  latter  two  reasons,  all  further  results  to  be  presented  will 
concentrate  on  the  stress  parallel  to  the  direction  of  heat  flow  for 
heating  rates  and  time  intervals  at  which  the  longitudinal  stresses 
exceed  the  transverse  stresses. 

Fig.  4  shows  the  distribution  of  the  thermal  stresses  along  the 

centerline  for  a  range  of  values  of  heating  rate  and  numerical  ve’ues 

-3  2 

for  the  thermal  diffusivity,  <  =  12.9  x  10  cm  /s,  width,  b  =  10  cm, 
length,  L  =  60  cm  at  the  instant  of  time,  t  =  1000  s.  The  magnitude  of 
stress  increases  linearly  with  heating  rate  with  the  position  of  peak 
stress  independent  of  heating  rate.  Nevertheless  with  increasing  heat¬ 
ing  rate  the  position  for  any  prescribed  value  of  stress  such  as  the 
tensile  fracture  strength  moves  towards  the  face  being  heated.  For  this 
reason,  the  location  of  tensile  failure,  if  it  occurs,  is  anticipated  to 
be  a  function  of  heating  rate. 

Fig.  5  for  a  value  of  heating  rate,  a  =  300°C/hr,  thermal  diffu- 
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sivity,  k  =  12.9  x  10  cm  /s,  and  width,  b  =  10  cm,  shows  the  stress 
distribution  along  the  centerline  for  a  range  of  values  of  time.  As 
expected,  the  magnitude  of  peak  stress  increases  with  increasing  time. 
Furthermore,  as  time  progresses  the  position  of  peak  stress  moves  away 
from  the  face  being  heated.  These  results  indicate  that  the  position  of 
failure  at  a  given  value  of  stress  is  expected  to  be  a  function  of  time. 
As  an  aside,  it  should  be  noted  that  at  the  highest  value  of  time 


indicated  in  Fig.  5  the  front  face  temperature  (417°C)  for  this  particu¬ 
lar  value  of  heating  rate  has  reached  a  value  at  which  in  many  candidate 
materials,  such  as  tar-bonded  magnesite  refractories,  substantial  soften¬ 
ing  and  creep  may  occur.  If  so,  considerable  relaxation  of  the  thermal 
stresses  may  occur  which  will  reduce  the  magnitude  of  the  thermal 
stresses  and  the  probability  of  failure. 

Fig.  6a  indicates  the  effect  of  the  magnitude  of  the  thermal  diffu- 
sivity  on  the  peak  value  of  longitudinal  tensile  thermal  stress  for 
three  values  of  time  and  a  heating  rate  a  =  300° C/hr,  width,  10cm, 
length,  L  =  60  cm.  These  data  indicate  that  the  magnitude  of  stress 
decreases  with  increasing  thermal  diffusivity.  This  effect  arises 
because  with  increasing  thermal  diffusivity  at  a  given  time  the  tem¬ 
perature  field  has  penetrated  to  a  greater  depth.  This  gives  rise  to  a 
greater  degree  of  temperature  uniformity  and  associated  lower  magnitude 
of  thermal  stress.  This  suggests  that  for  the  same  boundary  conditions 
in  practice  as  chosen  for  this  study  (linear  heating  of  the  front  sur¬ 
face)  ,  candidate  materials  should  have  a  thermal  diffusivity  as  high  as 
possible. 

Fig.  6b  shows  the  effect  of  the  value  of  thermal  diffusivity  for  t 
=  500  sec  in  greater  detail  for  the  values  of  heating  rate,  length  and 
width  identical  to  those  for  Fig.  6a.  These  results  indicate  that  the 
stresses  exhibit  a  maximum  at  an  intermediate  value  of  thermal  diffu¬ 
sivity,  also  barely  perceptible  for  t  =  500  sec  in  Fig.  6a.  This 
result  should  not  be  interpreted  as  suggesting  that  the  magnitude  of  the 
stresses  also  can  be  reduced  by  choosing  values  of  thermal  diffusivity 
as  low  as  possible.  This  would  be  true  only  for  a  value  of  thermal 


diffusivity  approaching  zero,  which  is  not  realizable  in  practice.  The 
decrease  in  stress  with  decreasing  thermal  diffusivity  at  the  lower 
values  of  thermal  diffusivity  is  because  at  the  value  of  times  for  these 
data,  the  temperature  field  has  not  been  able  to  penetrate  deeply  into 
the  refractory.  The  stresses,  however,  will  increase  in  magnitude  with 
increasing  time  as  clearly  indicated  in  Fig.  6a,  to  values  in  excess  of 
those  found  with  the  higher  values  of  thermal  diffusivity. 

Fig.  7  shows  the  effect  of  width  on  the  peak  value  of  tensile 
thermal  stress  for  a  range  of  values  of  thermal  diffusivity,  length  L  = 
60  cm,  t  =  500  s  and  heating  rate,  a  =  300°C/hr.  At  this  particular 
value  of  time  the  peak  value  of  stress  is  almost  independent  of  the 
thermal  diffusivity  as  shown  earlier  in  Fig.  6a.  The  stresses  exhibit 
a  maximum  at  an  intermediate  value  of  width.  This  latter  result  could 
be  of  practical  interest  as  it  would  permit  a  reduction  of  stress  in  a 
material  with  a  given  value  of  thermal  diffusivity  by  the  simple  expe¬ 
dient  of  increasing  or  decreasing  the  width,  depending  on  other  per¬ 
formance  criteria.  Additional  calculations  showed  that  the  magnitude  of 
the  peak  tensile  stress  is  independent  of  the  length  for  40  <  L  <  80  cm 
for  the  range  of  all  values  of  time,  width,  and  heating  rate  selected 
for  this  study.  For  this  reason,  it  does  not  appear  that  design  im¬ 
provements  to  reduce  the  probability  of  thermal  stress  fracture  can  be 
made  by  changes  in  length. 

All  the  results  presented  in  Figs.  2  to  7  were  based  on  the  assump¬ 
tion  of  plane  stress.  This,  in  effect,  treats  the  rectangular  structure 
as  a  thin  plate  in  which  the  stresses  in  the  direction  of  the  height  are 
identically  equal  to  zero.  In  practice,  however,  the  dimensions  of  the 


width  and  height  are  comparable  so  that  a  three-dimensional  effect  on 
the  stress  is  anticipated. 

A  three-dimensional  finite  element  model  consisted  of  rectangular 
solid  elements  with  five  elements  across  the  width,  four  elements  across 
the  height  and  twenty-five  elements  along  the  length.  As  in  the  plane 
stress  model  the  cross  section  element  dimensions  were  constant  along 
the  length,  but  the  element  length  dimension  was  much  smaller  at  the 
heated  face  gradually  getting  larger  toward  the  cold  face.  This  gives 
better  accuracy  where  the  stress  gradients  are  highest  near  the  heated 
face.  The  three-dimensional  analysis  was  limited  to  evaluation  of  the 
effect  of  cross  section  height  on  a  segment  for  the  same  numerical 

values  of  other  parameters  used  for  the  plane-stress  calculations: 
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diffusivity,  <  =  12.9  x  10  cm  / s,  width,  b  =  10  cm,  length,  L  =  60  cm 
and  heating  rate,  a  =  300°C/hr  at  a  time  of  5000  s.  The  values  of 
height  ranged  from  h  =  0  cm  (plane  stress)  to  10  cm,  which  corresponds 
to  a  square  cross  section.  Figure  8  shows  the  distribution  of  the 
stresses  along  the  centerline.  For  any  value  of  height  the  relative 
distribution  of  the  stresses  is  similar  to  the  plane-stress  distri¬ 
bution.  However,  the  peak  stress  magnitude,  increases  significantly 
with  increasing  height.  In  fact,  for  the  height  equal  to  the  width  of 
10  cm  the  stress  is  more  than  double  the  corresponding  value  for  plane- 
stress  conditions.  Fig.  9  shows  the  values  of  the  peak  stresses  plotted 
as  a  function  of  height.  The  shape  of  the  curve  suggests  the  plane 
stress  results  are  reasonably  valid  up  to  a  height  of  about  2.5  cm. 
Beyond  this  value  the  magnitude  of  stress  increases  almost  linearly  with 
height.  Of  course,  for  any  specific  geometry  encountered  in  practice 


involving  different  values  of  width  and  height  and  the  other  relevant 
parameters,  the  effects  will  need  to  be  established  independently.  In 
general,  however,  it  appears  that  geometries  approaching  square  cross 
sections  should  be  avoided.  In  fact,  it  appears  advisable  that  the 
value  of  height  of  some  5  cm  currently  encountered  in  practice  should  be 
reduced  coupled  with  an  increased  width  as  pointed  out  before. 

As  a  further  comment,  it  should  be  noted  that  linear  elastic 
stress-strain  behavior  was  assumed  throughout  for  all  data  presented 
above.  However,  varying  Young's  modulus  effects  are  expected  to  be 
significant,  especially  in  carbon- impregnated  or  tar-bonded  magnesia 
refractories.  Above  about  100°C,  Young's  modulus  undergoes  a  discon- 
tinous  change  from  about  70  to  11  GPa  as  the  direct  result  of  the  soften¬ 
ing  of  the  carbon  phase.  This  effect  on  the  numerical  results  of  this 
study  can  be  accounted  for  reducing  the  values  of  stress  by  the  same 
ratio  of  the  relative  reduction  of  Young's  modulus  if  the  whole  struc¬ 
ture  is  above  100°C.  In  principle,  this  effect  can  be  incorporated  in 
the  calculation  by  the  finite  element  method  regardless  of  the  magnitude 
and  distribution  of  temperature.  Such  calculations  were  judged  beyond 
the  scope  of  this  parametric  study,  but  should  be  carried  out  for  any 
design  of  given  geometry,  dimensions  and  other  pertinent  variables. 

Finally,  the  results  obtained  by  the  present  numerical  approach 
should  be  compared  with  those  predicted  by  the  theory  of  Kienow  ,  which 
states  that  the  stress,  a  along  the  centerline  parallel  to  the  length  is 
given  by: 

a  =  aE  (d2T/dx2)b2/(16  +  3b3/x3)  (5) 
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2  2 

where  d  T/dx  is  the  second  derivative  of  the  temperature  distribution 
at  the  position  x  behind  the  face  being  heated  and  all  other  symbols 
having  been  defined  earlier. 

Figs.  10a  and  10b  compare  the  distributions  of  stresses  along  the 

centerline  calculated  by  the  Kienow  theory  and  the  present  method  for  a 

-3  2 

value  of  thermal  diffusivity,  <  =  12.9  x  10  cm  /s,  two  different 
widths,  b  =  7.5  and  15  cm,  length,  L  =  60  cm,  at  the  instant  of  time,  t 
=  1000  s,  for  a  heating  rate  of  300°C/hr.  For  both  sets  of  data  the 
magnitude  of  the  peak  stresses  appear  to  be  in  reasonable  agreement.  At 
positions  other  than  those  for  the  peak  stress,  the  agreement  is  less 
good.  Such  disagreement  is  expected  to  be  particularly  severe  at  the 
rear  face  at  which  the  Kienow  theory  predicts  a  finite  value  of  stress, 
whereas  no  such  stress  perpendicular  to  a  free  surface  should  exist. 

When  the  effect  of  height  is  considered  the  Kienow  theory  underestimates 
stresses  by  a  large  factor.  For  this  reason,  it  is  advised  that  caution 
be  exercised  in  the  use  of  the  Kienow  theory  for  geometries,  dimensions, 
and  values  of  the  parameters  other  than  those  used  for  the  data  in  Fig. 
10. 

The  existence  of  the  stresses  calculated  above,  relies  on  the  non¬ 
linear  temperature  distribution  during  transient  heat-up.  However,  such 
stresses  can  exist  also  under  conditions  of  steady-state  heat  flow  if 
the  thermal  conductivity  of  the  material  of  construction  is  strongly 
dependent  on  temperature.  Indeed  this  is  the  case  for  the  many  mate¬ 
rials  in  which  heat  conduction  occurs  primarily  by  phonon  transport. 

The  magnitude  of  these  steady-state  stresses  probably  will  be  less  than 
those  encountered  during  transient  heat-up  unless  the  rate  of  heat-up  is 
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For  a  sample  calculation  of  c  the  following  values  will  be  assumed 

Ktt  =  2xl06,  MN*nf3/2,  v  =  0. 3,  a  =  13xlCf6  °C_1,  E  =  6xl010Nnf2 
lie 

4  -1 

and  e  =  10  °Cm  .  Substitution  of  these  values  into  Eq.  7.  yields: 

ccr  -  0.01  m  (1  cm)  (8) 

The  formation  of  a  crack  of  this  radius  following  the  initiation  of 
2 

fracture  is  not  unlikely  especially  in  relatively  dense  refractory 
materials  with  the  higher  values  of  tensile  fracture  stress.  Any  crack 
larger  than  c  will  be  unstable  and  will  propagate.  This,  in  turn, 
increases  the  crack  size  even  further,  which  enhances  the  degree  of 
instability.  Although  Eq.  7  is  strictly  valid  for  a  crack  in  an  in¬ 
finite  solid,  it  is  anticipated  that  in  a  segmented  structure,  the  crack 
will  propagate  to  the  outer  boundaries. 

This  phenomenon  is  offered  to  explain  the  observed  total  disintegra¬ 
tion  of  such  segmented  refractory  structures,  which  can  involve  the 
simultaneous  fracture  of  a  large  number  of  neighboring  segments,  with 
disastrous  consequences  for  the  total  structure  as  a  whole.  Possibly, 
this  failure  mode,  at  least  in  principle,  may  be  suppressed  by  reducing 
the  tensile  fracture  stress  as  much  as  possible  and  by  promoting  failure 

at  as  many  fracture  initiation  sites  as  possible.  As  indicated  theor- 
2 

etically  both  these  approaches  tend  to  decrease  the  crack  size  which 
results  following  the  initiation  of  failure.  Clearly,  this  approach  if 
practical,  also  will  depend  on  other  performance  criteria  imposed  on  the 
structure. 

Even  under  conditions  of  steady-state  heat  flow  with  much  lower 
values  of  temperature  gradient  than  assumed,  crack  instability  described 
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by  Eq .  7  can  exist.  This  can  lead  to  fatigue  failure  by  the  sub-critical 
growth  of  cracks  by  stress-corrosion  or  diffusional  mechanisms  depending 
the  environmental  variables  and  magnitude  of  temperature.  It  is  thought 
this  failure  mechanism  may  be  responsible  for  the  long-term  failure  of 
glass-tank  refractories  as  the  result  of  crack  formation  during  initial 
heat-up  as  described  in  this  paper. 

In  summary,  this  study  has  established  the  origin  of  the  mode  of 
thermal  stress  failure  of  segmented  structures  subject  to  rapid  heat-up. 
The  definition  of  the  role  of  the  various  pertinent  material  and  geo¬ 
metric  variables  hopefully  will  be  of  help  in  the  improved  design  of 
such  structures  in  engineering  practice.  Specifically,  the  results 
indicate  that  improvements  in  the  resistance  to  thermal  stress  failure 
of  segmented  refractory  structures  can  be  brought  about  by  increasing 
and  decreasing  the  values  of  the  width  and  height,  respectively,  from 
the  values  now  commonly  encountered  in  practice. 
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TABLE  1 


Properties  of  Magnesia-Based  Refractory  and  Parameter  Values 


Young's  Modulus  of  elasticity  (E): 

70.0  GPa 

Coefficient  of 
thermal  expansion  (a) 

12  x  10"6  0C_1 

Thermal  diffusivity  (k) : 

4.0  x  10"3  -  4  : 

Poisson's  ratio  (v) : 

0.3 

Length  (L) : 

40  -  80  cm 

Width  (b): 

5  -  18  cm 

Height  (h) : 

2.5  -  10  cm 

Heating  rate  (a): 

100  -  900  °C/hr 

Time  (t) : 


500  -  15000  s 


LIST  OF  FIGURES 


1  Schematic  illustration  of  the  three-dimensional  solid  segment  being 
heated  on  the  x  =  0  side  and  insulated  on  all  other  sides. 

2  Contour  plots  of  (a)  temperature,  (b)  longitudinal  stress,  and  (c) 

transverse  stress  in  a  segment  60  cm  long  by  15  cm  wide  for  a  heating 

-3  2 

rate,  a  =  300°C/hr,  a  thermal  diffusivity,  k  =  12.9  x  10  cm  /s, 
and  other  properties  in  Table  1  at  a  time  t  =  3000  s. 

3  Contour  plots  of  (a)  temperature,  (b)  longitudinal  stress,  and  (c) 

transverse  stress  in  a  segment  60  cm  long  by  15  cm  wide  for  a  heating 

-3  2 

rate,  a  =  300°C/hr,  a  thermal  diffusivity,  k  =  12.9  x  10  cm  /s,  and 
other  properties  in  Table  1  at  a  time,  t  =  500  s. 

4  Longitudinal  stress  distribution  along  the  centerline  for  a  range 

_3 

of  heating  rates  using  values  of  thermal  diffusivity,  k  *=  12.9  x  10 
2 

cm  /s,  width,  b  =  10  cm,  length,  L  =  60  cm  at  a  time,  t  =  1000  s. 

5  Longitudinal  stress  distribution  along  the  centerline  for  a  range  of 

-3  2 

times  using  values  of  thermal  diffusivity,  <  *  12.9  x  10  cm  /s,  heat¬ 
ing  rate,  a  =  300°C/hr,  width,  b  =  10  cm,  length,  L  =  60  cm. 

6  Peak  longitudinal  stress  as  a  function  of  thermal  diffusivity  and 
values  of  heating  rate,  a  =  300°C/hr,  width,  b  *  10  cm,  length, 

L  =  60  cm:  (a)  range  of  values  of  time,  and  (b)  expanded  scale 
for  t  =  500  s. 

7  Peak  longitudinal  stress  as  a  function  of  segment  width  for  three 
values  of  thermal  diffusivity  with  a  length,  L  =  60  cm,  a  heating 
rate,  a  =  300°C/hr  at  a  time,  t  =  500  s. 
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8  Longitudinal  stress  distribution  along  the  centerline  from  three- 

dimensional  analysis  for  a  range  of  heights,  d,  from  plane  stress 

to  a  square  cross  section  with  values  of  thermal  diffusivity, 

-3  2 

k  =  12.9  x  10  cm  /s,  length,  L  =  60  cm,  width,  b  =  10  cm,  heating 
rate,  a  =  300°C/hr  at  a  time,  t  =  5000  s. 

9  Peak  longitudinal  stress  as  a  function  of  height,  d,  from  plane 

stress  to  a  square  cross  section  with  values  of  thermal  diffusivity, 

-3  2 

k  =  12.9  x  10  cm  /s,  length,  L  =  60  cm,  width,  b  =  10  cm,  heating 
rate,  a  =  300°C/hr,  at  a  time,  t  =  5000  s. 

10  Longitudinal  stress  distributions  calculated  by  the  Kienow  equation 

and  finite  element  methods  for  a  length,  L  =  60  cm,  a  thermal  diffusivity, 
-3  2 

k  =  12.9  x  10  cm  Is,  a  heating  rate,  a=  "300°C/hr  at  a  time,  t  =*  1000  s, 
and  segment  widths:  (a)  b  =  15  cm,  (b)  b  =  7.5  cm. 
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Fig.  2.  Contour  plots  of  (a)  temperature,  (b)  longitudinal 
stress,  and  (c)  transverse  stress  in  a  segment  60  cm  long  by 
15  cm  wide  for  a  heating  rate,  a  ■  300°C/hr,  a  thermal 
diffusivity,  k  ■  12.9  x  10“  ^  cm^/s,  and  other  properties  in 
Table  1  at  a  time  t  ■  3000  s. 
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Fig.  3.  Contour  plots  of  Ca)  temperature,  (b)  longitudinal 
stress,  and  (c)  transverse  3tress  in  a  segment  60  cm  long  by 
IS  cm  wide  for  a  heating  rate,  a  *  300°C/hr,  a  thermal 
diffusivity,  <  ■  12.9  x  10"^  cm^/s,  and  other  properties  in 
Table  1  at  a  time,  t  ■  500  s. 


DISTANCE  FROM  HOT  FACE  ,  cm 


Fig-  4.  Longitudinal  stress  distribution  along  the 
centerline  for  a  range  of  heating  rates  using  values 
of  thermal  diffusivity,  k  =  12-9  x  10-^  cm^/s,  width, 
b  =  10  cm,  length,  L  *  60  cm  at  a  time,  t  =  1000  s. 


DISTANCE  FROM  HOT  FACE ,  cm 


Fig.  5.  Longitudinal  stress  distribution  along  the 
centerline  for  a  range  of  times  using  values  of 
thermal  diffusivity,  k  ■  12.9  x  10“^  cra^/s,  heating 
rate,  a  *  300°C/hr,  width,  b  =  10  cm,  length,  L  *  60  cm 
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LONGITUDINAL  STRESS,  MPa  PEAK  LONGITUDINAL  STRESS,  MPa 


(PLANE  STRESS ) 


HEIGHT,  cm 


Fig.  9.  Peak  longitudinal  stress  as  a  function  of 
height,  d,  from  plane  stress  to  a  square  cross  section 
with  values  of  thermal  diffusivity,  k  ■  12.9  x  10" 
cm  /s,  length,  L  =  60  cm,  width,  b  *  10  cm,  heating 
rate,  a  =  300°C/hr,  at  a  time,  t  =  5000  s. 


DISTANCE  FROM  HOT  FACE  ,  cm 


DISTANCE  FROM  HOT  FACE  ,  cm 


Fig.  10.  Longitudinal  stress  distributions  calculated 
by  the  Kienow  equation  and  finite  element  methods  for 
a  length,  L  =  60  cm,  a  thermal  diffusivity,  k  =  12.9  x 
10-3  cm3/ s ,  a  heating  rate,a=  300°C/hr  at  a  time, 
t  =  1000  s,  and  segment  widths:  (a)  b  e  15  cm,  (b) 

l_  _  *»  e 


7.5  cm. 
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ABSTRACT 

A  reviev  is  presented  of  recent  experimental  data  on  the  effect  of  micro- 
structural  and  compositional  heterogeneity  on  the  thermal  diffusivity/conducti- 
vity  of  structural  materials  for  high-temperature  applications.  These  materials 
include  nitride  ,  arb’aes,  oxides,  and  various  composites.  The  microstructural 
and  compos:',  otial  variables  include  microcracking,  solid  solution  alloying, 
densif ication  aios,  impurities,  and  the  amount,  type,  distribution,  and  orienta¬ 
tion  of  dispersed  phases.  Such  variables  can  be  used  to  tailor  the  heat  transfer 
properties  of  a  material  for  specific  requirements. 

INTRODUCTION 

The  increased  performance  requirements  for  many  components  of  aerospace 
structures,  energy-conversion  systems,  industrial  process  vessels,  etc.,  has  led 
to  an  increased  demand  for  materials  which  exhibit  superior  mechanical  and 
chemical  behavior  at  high  temperatures.  Such  materials  frequently  include  the 
refractor)'  nitrides,  carbides,  and  oxides. 

The  in-service  performance  of  such  materials  frequently  also  depends 
critically  on  their  thermal  properties  such  as  thermal  expansion,  and  thermal 
conductivity  and  diffusivitv.  Generally,  values  for  the  coefficient  of  thermal 
expansion  should  be  as  low  as  possible.  Likewise,  low  values  for  the  thermal 
conductivity  and  diffusivitv  are  required  for  good  thermal  insulation:  however, 
for  heat  exchangers,  the  thermal  conductivity  and  diffusivitv  should  be  as  high 
as  possible.  The  latter  requirement  also  is  critical  for  those  components  which 
are  subject  to  failure  by  thermal  stresses,  which  are  inevitably  encountered 
during  transient  and  steady-state  non-isothermal  conditions.  Thermal  stress 
failure  is  particularly  critical  for  manv  refractor)'  materials,  which,  because  of 
their  brittle  nature,  usually  undergo  totally  catastrophic  failure,  rendering  the 
component  incapable  of  continued  performance.  Frequently,  the  service  conditions 
for  such  materials  may  simultaneously  require  high  thermal  insulating  abilitv  and 
resistance  to  thermal  stress  failure.  In  such  cases,  appropriate  trade-offs  in 
materials  properties  need  to  be  made,  or  alternatively,  design  modifications  may 
be  required. 

Clearly,  for  the  design  and  selection  of  materials  for  high-tenperature 
service  conditions,  quantitative  information  for  the  values  of  the  thermal  con¬ 
ductivity  and  thermal  diffusivitv  of  candidate  materials  is  required,  together 


with  a  detailed  understanding  of  the  intrinsic  variables  which  may  affect  their 
magnitude. 

Heat  transport  through  solids  occurs  primarily  by  phonon,  photon  or  electron 
transport.  The  heat  flux  which  results  from  these  mechanisms  depends  strongly 
on  the  associated  specific  heat,  the  temperature,  and  the  existence  of  structural 
and  chemical  imperfections  such  as  vacancies,  dislocations,  grain  boundaries, 
foreign  atoms,  optical  discontinuities  and  other  variables  which  contribute  to 
phonon,  photon,  and  electron  scattering  (1,2) .  At  the  microstructural  level, 
the  conduction  of  heat  is  affected  by  the  presence  of  pores,  second-phase  inclu¬ 
sions,  cracks,  grain  boundary  phases,  texture,  and  other  factors. 

The  purpose  of  this  paper  is  to  present  an  overview  of  experimental  data 
obtained  by  the  present  writers  and  co-workers  over  the  last  few  years  which 
illustrates  many  of  the  above  microstructural  effects  on  heat  conduction.  The 
laser-flash  diffusivity  method  was  used  for  all  studies  ( 2)  ■  By  multiplying  the 
experimental  values  of  the  thermal  diffusivity  by  the  specific  heat  per  unit 
volume,  the  corresponding  values  for  the  thermal  conductivity,  when  required, 
were  readily  calculated.  For  many  of  the  experimental  details  on  the  material 
preparation  and  characterization,  the  reader  is  referred  to  the  original  publica¬ 
tions.  For  convenience,  the  data  will  be  presented  in  terms  of  a  number  of 
phenomena  observed  for  a  given  material  such  as  silicon  nitride,  or  else  in  terms 
of  a  given  phenomenon,  such  as  the  effect  of  microcracking,  observed  for  a  number 
of  different  materials. 

EXPERIMENTAL  DATA 

Silicon  Nitride 

Because  of  its  relatively  low  coefficient  of  thermal  expansion  and  other 
favorable  properties,  polycrystalline  silicon  nitride  is  a  prime  candidate  mate¬ 
rial  for  such  purposes  as  an  all-ceramic  turbine  engine  and  components  of  an 
adiabatic  diesel  engine.  Silicon  nitride  can  be  made  by  hot-pressing  appropriate 
powders  or  by  reaction  sintering,  which  consists  of  nitTiding  silicon  powders 
pressed  into  the  appropriate  shape. 

Hot-pressed  silicon  nitride  components  are  made  using  powders  which  consist 
initially  of  the  a-crystallographic  modification.  To  promote  densification,  a 
hot-pressing  aid  consisting  of  magnesium  oxide,  yttrium  oxide  or  other  oxide  is 
added  to  the  silicon  nitride  powder.  By  reacting  with  the  silicon  oxide  layer 
on  the  silicon  nitride  particles,  the  hot-pressing  aid  exists  in  the  form  of  a 
glassy  silicate  grain  boundary  phase.  At  the  hot  pressing  temperature,  the  glassy 
phase  promotes  densification  by  viscous  flow.  The  amount  of  hot-pressing  aid  can 
be  varied  in  order  to  optimize  the  properties  required  for  a  given  application. 
Since  the  sintering  aid  exists  as  a  separate  grain  boundary  phase,  the  heat  con¬ 
duction  behavior  of  hot-pressed  silicon  nitride  is  a  function  of  the  relative 
amount  of  hot-pressing  aid  present. 

Two  additional  effects  occur  during  the  hot-pressing  of  silicon  nitride. 
First,  the  <*-Si3N4  transforms  into  the  6-phase  by  a  time-dependent  dissolution- 
precipitation  process  via  the  viscous  grain  boundary  phase.  Secondly,  the  3- 
particles  become  rod-shaped  due  to  a  preferred  growth  direction  and  become 
oriented  such  that  the  axis  of  the  rod  is  perpendicular  to  the  direction  of 
hot-pressing  (A, 5) .  This  latter  effect  introduces  a  texture  which  results  in  an 
anisotropy  of  those  properties  sensitive  to  grain  geometry  and/or  crystallographic 
orientation. 

Figures  1  and  2  present  evidence  for  these  chree  effects  on  the  conduction 
of  heat  through  hot-pressed  silicon  nitride  as  found  by  Ziegler  and  co-workers 
(6,7).  Figure  1  shows  the  thermal  diffusivity  of  hot-pressed  silicon  nitride 
with  5  wt .Z  MgO  as  a  function  of  hot-pressing  time,  perpendicular  and  parallel 
to  the  hot-ptessing  direction.  The  increare  in  thermal  diffusivity  with  hot- 
pressing  time  results  from  two  causes.  At  the  shorter  values  of  time,  the 
increase  in  density(i .  e.,  reduction  in  relative  pore  content)  has  the  major  effect. 
At  the  longer  hot-pressing  time,  the  additional  increases  in  the  thermal  diffusi¬ 
vity  result  from  the  o-to-8  transformation,  which  is  essentially  complete  at 
approximately  3  hours.  It  is  believed  that  due  to  a  zone-refining  effect  during 
its  growth,  the  3-phase  attains  a  higher  level  of  puritv  and  consequently  a 
higher  thermal  diffusivity  than  the  o-phase.  Because  of  the  preferred  alignment 
of  the  axis  of  the  rod-like  8-particles,  the  thermal  diffusivity  perpendicular 
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Fig.  1  Thermal  diffusivity  of  hot-pressed  silicon  nitride 
at  room  temperature  as  a  function  of  hot-pressing 
time  and  direction  (after  ref.  6). 
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Fig.  2  Thermal  diffusivity  of  hot-pressed  silicon 
nitride  at  room  temperature  as  a  function 
of  MgO  content  and  hot-pressing  direction 
(after  ref.  6). 


Co  Che  hoc-pressing  direction  is  higher  chan  parallel  co  Che  hoc-pressing 
direction. 

Figure  2  shows  Che  data  for  the  thermal  diffusivicy  as  a  function  of  MgO 
content  of  silicon  nitride,  hoC-pressed  for  sufficient  duration  for  che  a-B 
transformation  to  be  complete.  At  approximately  331  MgO  the  thermal  diffusivity 
exhibits  a  maximum.  Below  32  MgO,  the  increase  in  thermal  diffusivity  can  be 
attributed  primarily  to  the  enhancement  of  density.  Above  32  MgO,  however,  the 
increasing  amount  of  glassy  grain  boundary  phase,  with  low  inherent  thermal 
conductivity,  leads  to  a  decrease  in  thermal  diffusivity  with  increasing  MgO  con¬ 
tent. 

A  significant  degTee  of  anisotropy  in  heat  conduction  behavior  of  silicon 
nitride  can  be  introduced  by  the  addition  of  second  phase  particles  such  as 
boron  nitride.  Such  S13N4-BN  composites  have  good  resistance  to  catastrophic 
fracture  by  thermal  shock.  Due  co  its  graphite-like  hexagonal  crystal  structure, 
boron  nitride  is  highly  anisotropic  in  many  of  its  properties.  Furthermore, 
boron  nitride  powders  exist  in  the  shape  of  flat  particles  with  the  long  dimen¬ 
sion  along  the  basal  plane.  During  the  powder  orientation  stage  during  the 
Initial  phases  of  hot-pressing,  che  BN  particles  take  on  a  preferred  orientation 
with  the  long  dimension  (i.e., the  basal  plane)  perpendicular  to  the  hot-pressing 
direction.  Boron  nitride  exhibits  the  highest  value  of  thermal  conductivity 
within  the  basal  plane.  Thus,  their  preferred  orientation  within  the  silicon 
nitride  leads  to  anisotropic  thermal  conductivity  of  the  hot-pressed  Si3Ni-BN 
composites.  Evidence  for  this  effect  was  found  by  Niihara  et,  al. (8)  and  is  shown 
in  Figure  3.  X-ray  analysis  revealed  that  che  presence  of  the  boron  nitride 
reduces  che  degree  of  orientation  of  the  S-S13N4  grains  which  diminishes  the 
anisotropy  effect  shown  in  Figure  1.  For  this  reason,  the  pronounced  anisotropy 
at  the  higher  values  of  BN  content,  as  indicated  by  the  data  in  Figure  3,  is 
primarily  due  to  the  preferred  orientation  of  the  boron  nitride  particles. 

Solid-solution  alloying  is  known  to  be  very  effective  in  reducing  the  phonon 
mean  free  path  in  dielectric  materials,  thus  producing  a  significant  decrease  in 
the  thermal  diffusivity  and  conductivity.  This  effect  is  demonstrated  by 
experimental  data  of  Lange  et.  al.  (9.)  in  Figure  4,  which  shows  the  effect  of 
alumina  alloying  on  the  thermal  diffusivity  of  silicon  nitride.  As  expected,  the 
relative  effect  is  larger  at  che  lower  levels  of  temperature  chan  at  the  higher 
temperatures. 

Silicon  Carbide 

SiC-based  materials  consisted  of  particulate  composites  of  SiC-Si  and  SiC- 
fiber  reinforced  glass-ceramics.  Composites  of  magnesium  oxide  and  silicon 
carbide  will  be  discussed  in  a  later  section  on  che  effect  of  microcracking  on 
thermal  conductivity. 

Silicon  carbide-silicon  composites  are  a  new  composite  material  for  high- 
temperature  turbine  applications.  Figure  5  shows  experimental  data  for  the 
thermal  diffusivity  of  such  SiC-Si  composites  at  room  temperature  as  a  function 
of  silicon  content  as  reported  by  Srinivasan  et.  aL  (10).  Over  the  range  of 
silicon  content  investigated,  the  thermal  diffusivity,  to  a  first  approximation, 
is  Independent  of  silicon  content.  It  is  noteworthy  that  the  thermal  diffusivity 
of  the  SiC-Si  composites  is  well  below  the  corresponding  value  for  the  silicon 
carbide.  From  a  practical  viewpoint,  these  SiC-Si  composites  are  much  better 
thermal  insulators  than  single-phase  silicon  carbide.  This  latter  observation 
is  of  vital  importance  for  turbine  applications  in  order  to  minimize  heat  losses 
and  thereby  improve  fuel  efficiency. 

Nevertheless,  the  much  lower  values  for  the  thermal  diffusivity  of  these 
composites  than  for  SiC  are  unexpected  since  reported  values  for  the  thermal 
diffusivity  of  silicon  metal  (U.)  are  comparable  if  not  larger  than  those  for 
silicon  carbide.  In  this  respect  then,  the  data  in  Figure  5  appear  not  to  be 
in  accordance  with  composite  theory  which  indicates  that  for  materials  with 
similar  values  for  specific  heat  and  density,  the  thermal  diffusivity  of  com¬ 
posites  should  be  intermediate  co  the  values  of  the  individual  components.  In 
order  to  resolve  this  apparent  discrepancy,  data  were  obtained  for  the  thermal 
diffusivicy  of  the  industrial  grade  silicon  metal  used  to  fabricate  these  com¬ 
posites  and  compared  with  the  data  for  a  high  purity  zone-refined  silicon.  The 
results  are  shown  in  Figure  6.  Clearly  the  industrial  grade  silicon  has  a  much 
lower  thermal  diffusivity  chan  the  high-purity  silicon.  It  is  interesting  to 
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Fig.  3  Thermal  diffusivity  of  silicon 
nitride-boron  nitride  composites 
at  room  temperature  as  a  function 
of  boron  nitride  content  and  hot- 
pressing  direction  (ref.  8). 
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Fig.  4  Effect  of  alumina  alloying  on  thermal 

diffusivity  of  silicon  nitride  (ref.  9). 
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Comparison  of  thermal  diffusivity  of  industrial  and 
zone-refined  silicon  (ref.  10). 


note  that  the  total  Impurity  content  in  the  industrial  silicon  was  no  more  than 
12,  with  boron  and  iron  as  principal  constituents.  Thus,  the  apparent  discrepancy 
of  the  composite  data  shown  in  Figure  5  is  explained,  and  the  critical  role  of 
impurities  in  governing  the  heat  conduction  behavior  of  these  types  of  materials 
is  clearly  evident. 

The  SiC-fiber  reinforced  glass-ceramics  consisted  of  a  crystallized 
lithium  aluminosilicate  glass  reinforced  with  fibers  of  an  amorphous  silicon 
carbide.  These  composites  exhibit  excellent  mechanical  behavior,  such  as  high 
fracture  toughness  at  elevated  temperatures.  Not  only  were  the  thermal  transport 
properties  of  these  composites  of  interest,  but  those  of  the  individual  silicon 
carbide  fibers  as  well.  At  this  time,  no  convenient  method  exists  by  which  the 
thermal  conductivity  of  non-conducting  dielectric  fibers  can  be  measured.  The 
approach  taken  to  assess  these  latter  values  was  to  calculate  the  fiber  conducti¬ 
vity  from  the  experimental  data  for  the  fibrous  composite  and  the  single-phase 
matrix  without  fibers,  by  means  of  the  theory  for  Che  thermal  conductivity  of 
composites  (12) . 

Figure  7  shows  the  experimental  data  obtained  by  Brennan  et.  al.  (13)  for  the 
thermal  diffusivity  of  the  glass-ceramic  matrix  without  fibers  and  the  composite 
with  49  vol.%  fibers,  oriented  parallel  and  perpendicular  to  the  direction  of 
heat  flow.  The  data  for  the  composites  do  not  differ  appreciably  from  those  of 
the  matrix,  which  indicates  that  the  values  for  the  thermal  conductivity  of  the 
glass-ceramic  matrix  and  the  fibers  are  comparable.  Figure  8  shows  the  values 
for  the  thermal  conductivity  and  diffusivity  of  the  fibers  calculated  from  the 
data  given  in  Figure  7,  and  independently  obtained  values  for  the  specific  heat 
and  density  of  the  matrix  and  fibers. 

Two  noteworthy  comments  can  be  made.  First,  the  fiber  conductivity  is 
comparable  to  that  of  amorphous  glasses,  with  the  value  of  thermal  diffusivity 
some  two  orders  of  magnitude  below  the  value  for  the  thermal  diffusivity  of 
crystalline  silicon  carbide  shown  in  Figure  5.  Secondly,  the  thermal  conductivity 
shows  a  positive  temperature  dependence,  in  contrast  to  the  strongly  negative 
temperature  dependence  exhibited  by  dielectric  solids,  especially  near  room 
temperature.  Both  these  effects  are  related  to  the  amorphous  nature  of  the 
silicon  carbide  fibers,  which  also  contain  appreciable  amounts  of  oxygen.  The 
amorphous  structure  of  the  fibers  limits  the  phonon  mean  free  path  to  exceedingly 
small  values,  which  results  in  very  low  thermal  conductivity  and  thermal  diffusi¬ 
vity.  These  values  are  reduced  to  such  an  extent  that,  as  suggested  by 
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Fig.  7  Thermal  diffusivity  of  lithium  aluaino-silicate 
(LAS)  glass-ceramic  and  LAS-sllicon  carbide 
fiber  composites  parallel  and  perpendicular  to 
fiber  direction  (ref.  13). 


the  positive  temperature  dependence  of  the  thermal  conductivity,  radiation  can 
make  a  significant  contribution  to  the  total  conductivity.  This  latter  conclu¬ 
sion  implies  that  the  thermal  conductivity  of  these  composites  should  increase 
with  specimen  thickness.  In  fact,  this  was  found  to  be  so.  This  effect  should 
be  taken  into  account  in  the  design  and  performance  analysis  of  structures  or 
components  siade  from  these  types  of  composites. 

Mlcrocracked  Materials 

Brittle  materials  with  high  densities  of  microcracks  offer  unique  properties 
such  as  very  low  Young's  modulus  of  elasticity  (19,15),  high  fracture  energy  (15) 
coupled  with  a  stable  mode  of  crack  propagation,  high  strain  at  fracture  (16,17), 
ease  of  machinability ,  and  excellent  thermal  shock  resistance  (18,19,20) . 

Tor  randomly  oriented  equal  sized  penny-shaped  cracks,  the  thermal  conducti¬ 
vity  (K)  is  given  by  (21) : 

K  -KJ1  +  8Mb3/9]_1  (1) 

where  K<j  is  the  thermal  conductivity  of  the  crack-free  solid,  R  is  the  number  of 
cracks  per  unit  volume,  and  b  is  the  crack  radius.  Equation  1  indicates  that 
microcracking  will  lead  to  a  decrease  in  thermal  conductivity,  as  expected,  since 
in  the  absence  of  radiative  heat  transfer,  the  cracks  act  as  barriers  to  heat 
flow.  Because  microcracking  has  little  or  no  effect  on  the  specific  heat  per 
unit  volume,  equation  1  is  expected  to  describe  the  effect  of  microcracking  on 
thermal  diffusivity  as  well.  6 

In  polycrystalline  single-phase  materials  microcracks  can  arise  from  the 
anisotropy  in  thermal  expansion  of  the  individual  grains.  In  brittle  composites, 
microcracks  can  result  from  mismatches  in  the  coefficients  of  thermal  expansion 
of  the  individual  components.  Brittle  structural  materials  generally  are  made 
at  high  temperature.  On  cooling  to  room  temperature,  the  non-uniform  thermal 
expansion  of  the  grains  (or  components  in  the  case  of  composites)  will  result  in 
the  creation  of  internal  stresses.  If  they  are  of  sufficient  magnitude,  these 
stresses  can  lead  to  microcracking.  This,  however,  requires  precursors  such  as 


Fig.  8  Thermal  conductivity  and  thermal  diffusivity  of 
amorphous  silicon  carbide  fibers  parallel  and 
perpendicular  to  the  fiber  direction  (ref  13). 


pores  at  grain  boundaries,  triple  points,  or  areas  of  less  than  complete  adhe¬ 
sion  between  the  individual  components.  Fundamentally,  microcracks  will  form  if 
the  combination  of  the  magnitude  of  internal  stress  and  the  size  of  the  crack 
precursors  is  such  that  the  resulting  stress  intensity  factor  equals  or  exceeds 
the  critical  stress  intensity  factor.  In  polycrystalline  materials,  the  size 
of  residual  pores  which  act  as  microcrack  precursors  generally  exhibits  a  one- 
to-one  correspondence  with  the  size  of  the  grains.  For  this  reason,  it  is  ex¬ 
pected  that  the  effect  of  microcracking  on  thermal  conductivity  or  thermal 
diffusivity  will  show  a  pronounced  dependence  on  grain  size.  Reheating  micro- 
cracked  materials  back  to  the  manufacturing  temperature  will  be  accompanied  by 
crack  closure  and  should  lead  to  a  positive  temperature  dependence  of  the  thermal 
conductivity.  At  sufficiently  high  temperature,  and  if  time  permits,  crack 
healing  can  occur  by  diffusional  processes.  This  latter  effect  would  lead  to  a 
thermal  conductivity  which  i6  a  function  of  thermal  history.  Evidence  for  these 
effects  will  be  shown  in  the  subsequent  figures. 

Figure  9a  shows  data  measured  by  Siebeneck  et. al. (22)  for  the  temperature 
dependence  of  the  thermal  diffusivity  of  an  as  hot-pressed  polycrystalline  iron 
titanite  with  a  mean  grain  size  of  approximately  1  urn.  These  data  indicate  the 
normal  negative  temperature  dependence  of  a  dielectric  solid  in  which  phonon  con¬ 
duction  is  the  primary  mechanism  of  heat  transfer.  Figure  9b  shows  the  thermal 
diffusivity  of  an  iron  titanate  sample  which  was  given  an  annealing  treatment 
following  hot-pressing  in  order  to  promote  grain  growth.  The  resulting  grain 
size  distribution  was  bimodal  with  grains  as  large  as  20-30  pm  contained  within 
the  original  fine  grained  matrix.  Comparison  of  Figure  9a  and  9b  indicates  that 
the  thermal  diffusivity  of  the  large  grained  iron  titanate  has  only  about  one 
third  of  the  value  of  the  fine  grained  material.  This  effect  is  due  to  the 
microcracking  in  the  coarse-grained  material.  Figure  9b  also  indicates  the 
existence  of  a  pronounced  hysteresis  effect  as  the  specimen  is  heated  and  cooled. 
During  the  heating  part  of  the  cycle  the  increase  in  the  thermal  diffusivity  is 
attributed  to  crack  closure.  The  higher  values  obtained  during  the  cooling  part 
are  thought  to  be  due  to  either  irreversible  crack  closure  due  to  frictional 
effects  or  to  crack  healing  due  to  diffusional  effects.  Both  these  effects  will 
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Fig.  9  Thermal  diffusivity  of  polycrystalline 
iron  titanate.  a:  grain  size  of  -  1 
urn;  b:  grain  size  up  to  30  um  (ref.  22). 


contribute  to  the  observed  hysteresis  since  a  finite  degree  of  cooling  is  required 
before  healed  microcracks  will  repropagate.  As  Indicated  by  the  discontinuity  in 
slope  in  the  cooling  part  of  the  curve,  this  recracking  occurs  at  a  temperature 
of  approximately  400  C.  Such  crack  healing  and  re-formation  quite  likely  is 
accompanied  by  permanent  microstructure  changes,  and  the  crack  geometry  and 
dimensions  are  not  likely  duplicated  exactly  from  cycle  to  cycle.  Regardless  of 
whether  this  latter  effect  occurs  or  not,  Figure  9b  clearly  indicates  that  the 
effect  of  microcracking  on  thermal  conductivity  or  thermal  diffusivity  depends  on 
prior  thermal  history. 

Additional  aspects  related  to  the  effect  of  microcracking  on  thermal  diffusi¬ 
vity  are  shown  in  Figs.  10a  and  10b  for  a  fine  and  coarse  grained  magnesium  di- 
titanate,  respectively,  studied  by  Siebeneck  et.  aL  (23).  In  Figure  10a,  the  data 
for  the  thermal  diffusivity  during  cooling  are  well  below  those  obtained  during 
heating.  It  is  thought  that  this  effect  arises  from  the  diffusional  growth  of 
the  oicrocrack  precursors,  which,  during  the  initial  cooling  from  the  hot-pressing 
temperature,  were  not  of  sufficient  size  to  propagate  into  microcracks.  The 
additional  time  at  the  higher  temperatures  during  the  measurement  of  the  thermal 
diffusivity  permitted  sufficient  growth  so  that  on  cooling,  microcracks  formed 
which  were  not  present  during  the  Initial  cooling. 
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Fig.  10  Thermal  diffusivity  of  a.  fine  grained 

and  b.  coarse  grained  magnesium  dltitanate 
(ref.  23). 


Figure  XOb  shows  evidence  of  the  combined  effect  of  microcrack  healing  and 
precursor  growth  at  the  higher  temperature  for  a  coarser-grained  (70  um)  magesium 
dititanate  (23),  which  results  in  a  reversal  in  the  direction  of  hysteresis  at 
the  lower  and  higher  tenperature  changes.  Again,  the  prior  thermal  history  of 
the  specimen  can  have  a  major  influence  on  the  effect  of  microcracking  on  heat 
conduction  behavior.  Clearly,  for  microcracked  materials,  no  handbook  values  for 
thermal  conductivity  or  diffusivity  can  be  established. 

Experimental  data  for  microcracked  composites  permit  a  quantitative  evalua¬ 
tion  of  the  effect  of  microcracking  on  Che  thermal  conductivity.  Figure  11  com¬ 
pares  che  measured  values  for  the  thermal  diffusivity  of  a  microcracked  composite 
which  consists  of  a  magnesium  oxide  matrix  with  a  dispersed  phase  of  silicon 
carbide  (24),  with  values  calculated  from  composite  theory,  which  implicitly 
assumes  that  microcracks  are  absent  (12) .  Due  to  the  presence  of  the  microcracks, 
the  experimental  data  fall  well  below  the  composite  data.  The  difference  between 
these  values  increases  with  increasing  SiC  concent.  This  is  expected  since  the 
microcracks  in  the  matrix  are  associated  with  the  silicon  carbide  inclusions. 

At  the  highest  values  of  SiC  content,  the  microcracks  have  decreased  the  thermal 
diffusivity  by  nearly  a  factor  of  three. 

The  MgO-SiC  composites  also  show  evidence  for  crack  healing  at  the  higher 
temperatures  and  re-formation  on  cooling  as  3hown  in  Figure  12  for  a  composite 
with  30  wt.%  SiC.  Crack  re-formation  appears  to  occur  over  the  temperature 
range  of  about  850  to  600  C. 

Other  Materials 

An  interesting  effect  is  shown  by  composites  of  a  sodium  borosilicate  glass 
with  Ni  inclusions  studied  by  Powell  and  co-workers  (25) .  The  composites  were 
fabricated  in  order  to  study  the  composite  behavior  resulting  from  lack  of  adhe¬ 
sion  between  the  glass  matrix  and  Ni  inclusions.  These  composites  were  made  by 
hot-pressing  mixtures  of  powders  of  the  appropriate  glass  and  spherical  Ni 
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Fig.  11  Comparison  of  measured  and 
calculated  thermal  iffusi- 
vity  at  room  temperature  of 
microcracked  composites  of 
magnesium  oxide  matrix  con¬ 
taining  silicon  carbide 
(ref.  24). 
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Fig.  12  Temperature  dependence  of  ther¬ 
mal  diffusivity  of  MgO-SiC 
composite  with  30%  SiC  (ref. 
24). 


particles  at  a  temperature  of  approximately  700  C.  Since  the  coefficient  of 
thermal  expansion  of  the  Ni  is  much  higher  than  that  of  the  glass,  the  NI 
particles  separate  from  the  glass  phase  during  cooling,  which  results  in  the 
formation  of  an  interfacial  gap.  Such  a  gap  prevents  direct  heat  flow  through 
the  Hi  particles,  preventing  them  from  contributing  to  the  heat  conduction  pro¬ 
cess,  with  the  exception  of  some  heat  flow  through  the  point9  of  contact  between 
the  Ni  particles  and  the  glass  matrix. 

Figure  13  compares  the  experimental  data  for  the  thermal  diffusivity  of  these 
composites  with  a  range  of  Ni  contents  with  those  calculated  from  composite 
theory  (^) .  The  data  clearly  indicate  that  near  room  temperature  the  interfacial 
gap  lowers  the  thermal  diffusivity  appreciably.  This  effect  is  largest  for  the 
composite  with  the  highest  Ni  content,  as  expected.  However,  evidence  of  the 
Curie  point  indicates  that  the  nickel  still  contributes  to  the  thermal  diffusivity 
in  spite  of  the  interfacial  gap.  It  is  thought  that  this  effect  arises  from  the 
points  of  contact  between  the  Ni  and  the  glass,  which  still  permit  heat  flow 
into  the  Ni.  In  this  way  the  Ni  contributes  to  the  specific  heat  but  not  to  the 
conduction  of  heat  in  this  temperature  range.  At  temperatures  from  about  300  C 
to  600  C  the  thermal  diffusivity  shows  a  rapid  increase  with  temperature,  attri¬ 
butable  to  gap  closure  and  increased  contact  between  the  glass  and  nickel.  Such 
a  large  positive  temperature  dependence  i3  rather  uncommon  in  most  materials  over 
this  temperature  range.  The  mechanism  involved,  at  least  in  principle,  permits 
the  temperature  dependence  of  the  thermal  conductivity  or  diffusivity  to  be 
tailon.J  for  specific  purposes,  if  desired. 

For  dt"ign  purposes,  a  thermal  conductivity  of  dielectric  materials  which  is 
independent  of  temperature  may  be  required.  For  instance,  this  is  the  case  in 
linear  heat  flow  which,  for  a  constant  thermal  conductivity,  yields  a  linear 
temperature  gradient  and  zero  thermal  stress.  This  aim  can  be  achieved  by  syn¬ 
thesizing  composites  in  which  the  negative  temperature  dependence  of  one  component 
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Fig.  13  Comparison  of  calculated  and  measured 
values  of  che  thermal  diffusivity  of 
borosillcate  glass  with  Ni  inclusions 
(ref.  25). 


is  compensated  by  Che  positive  temperature  dependence  of  the  other.  A  specific 
composite  'or  which  this  aim  was  achieved  consisted  of  a  continuous  sodium 
borosilicate  glass  matrix  with  a  dispersed  phase  of  crystalline  aluminum  oxide. 
Above  room  temperature  the  former  material  exhibits  a  thermal  conductivity  which 
increases  with  temperature,  whereas  for  the  aluminum  oxide  over  the  same  tempera¬ 
ture  range,  the  chermal  conductivity  decreases  rapidly  with  temperature. 

Figures  14a  and  14b  show  the  experimental  data  for  the  chermal  diffusivity 
and  the  calculated  values  for  the  thermal  conductivity  of  this  composite  for  a 
range  of  alumina  contents  as  determined  by  Hasselman  et.al.(26).  At  approximately 
35  wt.Z  AI2O3  the  thermal  conductivity  is  virtually  constant  from  room  tempera¬ 
ture  to  600  C. 

FINAL  REMARKS 

the  data  presented  in  this  paper  clearly  indicate  the  important  role  of 
composition  and  microstructure  in  governing  the  heat  transfer  characteristics  of 
brittle  structural  materials.  In  general,  this  will  require  that  for  each  new 
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Ftg.  1.4  Heat  conduction  behavior  of  sodium  borosilicate  glass- 
alumina  composites,  a.  Experimental  values  for  the 
thermal,  diffusivit  •  and  b.  data  for  the  thermal  con¬ 
ductivity  calculated  from  the  thermal  diffusivity  (ref.  16). 
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batch  of  materials  the  values  for  the  thermal  conductivity  and  diffuslvlty  be 
measured,  rather  than  relying  on  handbook  values.  On  the  other  hand  the  high 
sensitivity  of  the  thermal  conductivity  and  dif fusivity  to  a  wide  range  of 
compositional  and  raicrosu'uctural  variables  provides  the  materials  technologist 
with  a  powerful  tool  to  tailor  these  properties  to  meet  specific  design  require¬ 
ments. 
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